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Abstract 

The description of nuclear reactions induced by supernova neutrinos has witnessed significant 
progress during the recent years. On one hand this progress is due to experimental data which 
serve as important constraints to model calculations, on the other hand it is related to advances 
in nuclear modelling itself and in computer hardware. At the energies and momentum transfers 
relevant for supernova neutrinos neutrino-nucleus cross sections are dominated by allowed tran¬ 
sitions, however, often with non-negligible contributions from (first) forbidden transitions. For 
several nuclei allowed Gamow-Teller strength distributions could be derived from charge-exchange 
reactions and from inelastic electron scattering data. Importantly the diagonalization shell model 
has been proven to accurately describe these data and hence became the appropriate tool to calcu¬ 
late the allowed contributions to neutrino-nucleus cross sections for supernova neutrinos. Higher 
multipole contributions are usually calculated within the framework of the Quasiparticle Random 
Phase Approximation, which describes the total strength and the position of the giant resonances 
quite well. Both are the relevant quantities for a reliable calculation of the forbidden contributions 
to the cross sections. 

The current manuscript reviews the recent progress achieved in calculating supernova-relevant 
neutrino-nucleus cross sections and discusses its verification by data. Moreover, the review sum¬ 
marizes also the impact which neutrino-nucleus reactions have on the dynamics of supernovae 
and on the associated nucleosynthesis. With relevance to the supernova dynamics these include 
the absorption of neutrinos by nuclei (the inverse of nuclear electron capture which is the domi¬ 
nating weak-interaction process during collapse), inelastic neutrino-nucleus scattering and nuclear 
de-excitation by neutrino-pair emission. For supernova nucleosynthesis we discuss the role of 
neutrino-induced reactions for the recently discovered vp process, for the r-process and for the 
neutrino process, for which neutrino-nucleus reactions have the largest impact. Finally we briefly 
review neutrino-nucleus reactions important for the observation of supernova neutrinos by earth- 
bound detectors. 
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1 Introduction 


February 1987 was the birth of extrasolar neutrino astronomy when detectors in Japan and the United 
States registered neutrinos which travelled over 50 kpc from the Large Magellanic Cloud to Earth and 
gave the hrst indications that a star in this neighbor galaxy of the Milky Way had exploded as a 
supernova [1, 2], This extraordinary scientific event also proved the general expectation that neutrinos 
are produced in enormous numbers in supernovae triggered by the core collapse of massive stars. In 
fact, about 99% of the gravitational binding energy released in the cataclystic event is carried away by 
neutrinos, clearly overpowering the kinetic energy associated with the expansion of the supernova and 
the energy radiated away as light, despite the fascinating fact that supernovae can shine as bright as 
entire galaxies. 

Although the neutrinos observed from supernova SN1987A were likely all electron antineutrinos, 
identihed by the Cerenkov light produced by the relativistic positrons after a charged-current neutrino 
reaction on the protons in the water Cerenkov detectors, the amount of observed neutrinos and their 
energy spectrum (of order a few 10s of MeV) conhrmed the general understanding of supernova dynam¬ 
ics [3]. These observations were supplemented by detailed studies of the SN1987A lightcurve, which as 
expected, followed the sequence of half-lives of radioactive nuclides like ^®Ni, ®^Ni or ^“^Ti, which were 
copiously produced in the hot supernova environment [4]. 

Core-collapse, or Type II, supernovae are the final fate of massive stars when at the end of hydrostatic 
burning their inner core, composed of nuclei in the iron-nickel mass range, runs out of nuclear fuel 
and collapses under its own gravity triggering an explosion during which most of the star’s material, 
partly processed in the hot environment in what is called explosive nucleosynthesis, is ejected into 
the Interstellar Medium [3]. In the general picture of core-collapse supernovae and their associated 
nucleosynthesis, neutrino reactions on nucleons and nuclei play an important role. Arguably, however, 
the most important impact neutrinos have on the supernova dynamics is the fact that, for the typical 
supernova neutrino energy scales, they virtually do not interact with matter for densities smaller than a 
few 10^^ g cm“^. This makes neutrinos the most efficient cooling mechanism during the late hydrostatic 
burning stages and the early collapse phase, where neutrinos could be thought of as free streaming 
with an appropriate energy loss correction for nuclear reactions mediated by the weak interaction like f3 
decay or electron capture. This picture does not hold at higher densities, say in excess of 10^^ g cm“^, 
the neutrino interaction with matter is strong enough making a detailed bookkeeping of neutrinos and 
their interaction with matter a tedious, but necessary requirement for supernova modeling. During the 
collapse phase when densities larger than 10^^ g cm“^ are reached, it is the elastic scattering of neutrinos 
on nuclei which changes the neutrino transport through the dense matter to a diffusion problem with 
time scale larger than the competing collapse time scale. As a consequence neutrinos are trapped during 
the final stage of the collapse. By inelastic scattering on electrons and, to a lesser extent, on nuclei, 
neutrinos exchange energy with matter and get thermalized. Inelastic neutrino-nucleus scattering plays 
also an interesting role in a short episode after core bounce, where it alters the spectrum of electron 
neutrinos emitted during the so-called neutrino burst [5]. 

The supernova explosion is triggered by a shock wave which, when passing outwards through the 
Fe-Ni-core, dissociates the heavy nuclei into free nucleons. The interaction of neutrinos, produced by the 
hot matter of the freshly-born neutron star in the center, with the free protons and neutrons behind the 
shock are an effective addtional energy source which, together with effects like convection and plasma 
instabilities, are required for successful explosions, as modern multi-dimensional supernova simulations 
show [6, 7]. 

It is also the competition of the various interactions of electron neutrinos with neutrons and of anti¬ 
electron neutrinos with protons, which determines the proton-to-neutron ratio of the matter ejected 
from the surface of the nascent proto-neutron star crucially influencing the subsequent nucleosynthesis 
of the ejected matter in this neutrino-driven wind model. Simulations indicate that there are periods 
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where the neutrino-driven wind matter is proton-rich (with an electron-to-nucleon ratio Y^. > 0.5) 
and where it is neutron-rich (Y). < 0.5). In the former case the nucleosynthesis occuring when the 
ejected matter reaches cooler regions at larger distances from the neutron star surface gives rize to the 
recently discovered vp process. For several years the ejection of neutron-rich matter in the neutrino- 
driven wind has been considered as the favorite site for the astrophysical r-process. However, modern 
supernova simulations predict astrophysical conditions which allow for a ’weak r-process’ [ 8 ], which 
can produce elements up to the second peak around mass number A ~ 130, which are, however, not 
sufficient for the production of r-process nuclides at the gold-platinum peak and the transactinides. As 
the nucleosynthesis in the neutrino-driven wind scenario occurs in the presence of enormous neutrino 
fluxes, neutrino-induced reactions on nuclei might play an important role in these processes. Particularly 
relevant are here neutrino-induced reactions with particle emission (neutron, protons, a) in the hnal 
channel as they implicitly alter the nuclear abundance distributions. Neutrino-induced particle emission 
of nuclei is also the key process in the neutrino nucleosynthesis process [9] which is responsible for the 
production of selected nuclides generated by particle spallation of more abundant nuclides in charged- 
or neutral-current reactions in the outer shells of the star. 

Obviously neutrino-induced reactions by the various target nuclei are the means for detection of 
supernova neutrinos. After the observation of anti-electron neutrinos from SN1987A, it is the aim 
of operational and future detectors to individually detect the various neutrino types and proof the 
differences in the spectra of v^. and (usually used to combine v^^ and neutrinos and their 
anti-particles which are expected to have quite similar spectra) predicted by supernova models. 

Neutrino reactions on nuclei can be described in perturbation theory. For the energies and momen¬ 
tum transfers relevant for supernova neutrinos the task reduces to a nuclear structure problem where 
mainly allowed and hrst-forbidden transitions are of importance. Fermi transitions, of relevance for 
elastic neutrino scattering and for (i/g, e“) reactions, can only occur to Isobaric Analog States. Due 
to progress in computer hardware and impressive advances in shell model developments Gamow-Teller 
(GT) strength distributions can now be calculated for up to medium-mass nuclei taking the relevant 
nuclear degrees of freedom into account. This theoretical progress has been accompanied by novel 
experimental techniques which allow the determination of the GT_ strength (in which a neutron is 
changed into a proton as it is relevant for (z/e,e“) reactions) by (p, n) [10] and (^He,f) [11, 12] charge- 
exchange reactions and of the GT+ strength (in which a proton is changed into a neutron as in the 
(z/e,e+) reaction or in electron capture) by (n,p) [13] and, with significantly improved resolution, by 
the (d,^ He) [14] and (t,^He) [15] reactions at intermediate energies. The GTq strength, important for 
inelastic neutrino-nucleus scattering, can be studied in (p,p’) reactions, but it is also constrained by 
precision Ml data obtained in inelastic electron scattering for spherical nuclei where the orbital con¬ 
tribution to the Ml transition is small and the spin part is dominated up to a constant factor by the 
same operator as the GTq transitions. Importantly, the shell model studies reproduce the measured 
GT and Ml strength distributions quite well, except for a constant normalization factor and hence 
it is the method of choice to describe the Gamow-Teller contributions to the neutrino-nucleus cross 
sections. The calculation of first-forbidden (dipole) transitions require multi-shell studies which, due to 
computer limitations, have only been possible for light nuclei yet. Usually dipole contributions to the 
neutrino cross sections are derived within the Quasi-Particle Random Phase Approximation (QRPA) 
which can handle model spaces with a larger number of single particle orbitals than the shell model 
space, however, for the price of considering only (2 particle - 2 hole) correlations among nucleons. The 
QRPA is well suited to describe the energy centroid and total strength of giant resonances and hence 
gives a fair description of neutrino-nucleus cross sections if the neutrino energy is sufficiently large so 
that the cross sections are dominated by the contributions from giant resonances. 

Neutrino-induced spallation reactions are usually described in a two-step process, calculating the 
neutrino-induced excitation spectrum within the framework of the shell model and QRPA, followed by 
a determination of the decay probability into the various particle (or 7 ) channels using the statistical 
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model. Exploiting this two-step strategy differential cross sections have been derived for nearly the 
entire nuclear chart which can now been routinely used in nucleosynthesis studies. 

In this article we would like to review recent developments in the description of neutrino-nucleus 
reactions. Hereby we will focus solely on neutrino energies up to a few 10s of MeV, as they are relevant 
for supernovae. The description of reactions induced by higher-energetic neutrinos as they are important 
for atmospheric, accelerator or astrophysical neutrinos from other sources require different theoretical 
approaches; as outlined e.g. in Refs. [16, 17]. Our review can be viewed as an extension and update 
of the article by Kolbe et al. [18]. Hence we will only briefly repeat some of the material presented 
there, but focus mainly on the research performed later. Our review is organized as follows: in the next 
section we will briefly describe the current understanding of core-collapse supernovae, paying special 
attention to the role of neutrinos. In particular we will discuss the neutrino spectra during the collapse 
and arising from the cooling of the proto-neutron star. In Section 3 we summarize the basic cross section 
formulae for neutrino-nucleus reactions and discuss the models, mainly shell model and QRPA, used to 
derive the necessary nuclear input. We will also compare strength functions obtained with these models 
with experimental data. Furthermore we discuss in Section 3 on the basis of the statistical model the 
derivation of the decay probabilities into particle and 7 channels as required to calculate partial cross 
sections for neutrino-induced spallation reactions and summarize some general results. Sections 4 and 5 
are concerned with the description of specific neutrino-nucleus reactions and their role for the supernova 
dynamic and nucleosynthesis, respectively. Finally we summarize in Section 6 the progress achieved to 
describe neutrino-nucleus reactions of importance for operational and future supernova detectors. 


2 Brief description of core-collapse supernovae 

During their long lifetimes of millions and billions of years stars generate the energy necessary to achieve 
hydrostatic equilibrium by nuclear reactions in their interior [19, 20, 21]. The first nuclear reaction 
source which a star tabs is hydrogen burning, the fusion of 4 proton nuclei to ^He. Hydrogen burning 
occurs, in stars like our Sun, in a reaction sequence called pp-chains beginning with the formation of a 
deuteron by the fusion of two protons mediated by the weak interaction. In more massive stars than the 
Sun hydrogen burning occurs at slightly higher temperatures making hydrogen burning via the CNO 
cycle energetically more efficient than by the pp chains [22, 23]. It is interesting to note that hydrogen 
burning produces a nucleus (^He) with an identical amount of protons and neutrons (i.e. W = 0.5) 
from matter with an unbalanced proton-to-neutron ratio basically inherited from the Big Bang. It is 
an important feature of hydrostatic burning and nuclear stability that this value of W = 0.5 is kept 
during the various burning stages. 

The ashes of hydrogen burning, ^He, can be burnt once the temperature of the star’s interior has 
increased sufficiently by core contraction to enable the fusion of three “^He nuclei to This triple-alpha 
reaction and the subsequent ^^C(q;, 7)^®0 reaction constitute helium burning. Both together determine 
the carbon/oxygen ratio in the Universe. Hence helium burning is essential for life as we know it. 

In sufficiently massive stars the sequence of core contraction associated with a rise of the temperature 
repeats itself igniting successively the ashes of the previous burning phase. Thus following helium 
burning massive stars undergo in their centers the sequence of carbon, neon, oxygen and silicon burning 
until an ’iron core’ is produced as the product of silicon core burning. At this phase of its hydrostatic 
burning, a massive star consists of concentric shells that are the remnants of its previous burning phases 
(hydrogen, helium, carbon, neon, oxygen, silicon) surrounding the iron core. 

The temperature in the iron core (a few 10® K) is sufficiently high to establish an equilibrium 
of reactions mediated by the strong and electromagnetic interaction and their inverse. Under such 
conditions the matter composition is given by Nuclear Statistical Equilibrium (NSE). Importantly once 
such an equilibrium is achieved nuclear reactions (by strong and electromagnetic forces) cease as energy 
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sources. Hence the star has lost its nuclear energy source in the center and the iron core becomes 
eventually unstable under its own gravity. We note that the NSE matter composition depends on the 
astrophysical parameters (temperature and density). For the following discussion it is important that 
it also depends on the proton-to-neutron ratio, i.e. on the W value, which can only be changed by 
reactions mediated by the weak interaction. Initially the ashes of silicon burning consists roughly of 
the same amount of protons and neutrons (Yg ~ 0.5) favoring, in NSE, nuclei in the vicinity of the 
double-magic nucleus ^®Ni. Under (early) core conditions reactions mediated by the weak interaction 
are not in equilibrium with their inverse. This, as we will discuss next, will change the Yg value during 
the collapse and relatedly the matter composition. 

While the mass of the core is dominated by the nuclear component, the pressure is given by the elec¬ 
trons, which under core conditions, form a degenerate ultrarelativistic Fermi gas; both components are 
connected by the electron-to-nucleon ratio Yg. Once the core density reaches values of order 10® g cm“^ 
or higher, the electron Fermi energy grows to values in excess of MeV, the core matter can efficiently 
lower its free energy by electron captures on nuclei. This process has three important consequences: 

1) At densities below 10^^ g cm“^ the neutrinos produced by the capture process can leave the star 
unhindered carrying away energy. This is a very efficient cooling mechanism which, as noted in the 
fundamental paper by Bethe, Brown, Applegate and Lattimer (BBAL) [24] keeps the entropy of the 
matter low. As a consequence heavy nuclei survive during the collapse. 

2 ) Electron capture reduces the number of electrons (i.e. Yg gets smaller). As a result the pressure 
which the matter can stem against the gravitational collapse is reduced and the collapse is accelerated. 

3) As electron captures change a proton in the nucleus into a neutron, the matter composition is 
driven towards more neutron-rich (and heavier) nuclei. With increasing neutron excess, nuclei become 
unstable against [3 decay. Under core conditions these decays are, however, hindered due to Pauli 
blocking by electrons and at sufficiently high densities (say in excess of 10^® g cm“®) the phase space 
for beta decays closes and electron capture dominates by a large factor driving the core composition to 
neutron-rich, /^-unstable nuclei 

It should be mentioned that f3 decay is temporarily competitive with electron capture after silicon 
depletion and during silicon shell burning [25, 26]. Besides counteracting electron captures by increasing 
the Yg value, f3 decays are also an important neutrino source and thus add to the cooling of the core and 
the reduction of the entropy [25]. Finally, we remark that it is quite fortuitous that for pre-supernova 
simulations {3 decay rates are only required for nuclei for which large-scale shell model calculations can 
be performed [27]. 

Exploiting progress in nuclear modelling electron capture rates at supernova conditions have recently 
been calculated using an appropriately chosen hierarchy of approaches [28]. At densities up to about 
IQio g supernova weak-interaction rates are dominated by nuclei in the mass range A = 45 — 65 
(p/-shell nuclei). At such densities the electron chemical potential is of the same order as the nuclear 
Q values requiring a detailed reproduction of the GT+ distribution due to the strong phase space 
sensitivity of the rates. As we will discuss in the next Section, modern shell model calculations are 
capable to describe GT_|_ distributions rather well. Shell model diagonalization in the pf model space 
has therefore been used to calculate the weak-interaction rates for the nuclei with A ~ 50—65 [29, 27, 30]. 
These rates explicitly consider that, due to the hnite temperature in the star, electron capture can also 
occur from excited states, which can have different GT+ distributions than the ground state. With 
continuing electron captures and at higher densities (above 10^® g cm“®) nuclei become relevant for 
the capture rates which have proton numbers Z < 40 and neutron numbers N > 40; i.e. in the 
Independent Particle Model (IPM) the valence protons are in the pf shell, while the valence neutrons 
occupy already the next major shell {sdg shell). For such nuclei the IPM predicts a vanishing GT+ 
strength as the GT operator acts only in spin-isospin space and all GT+ transitions of the valence 
protons would be blocked due to the complete occupation of pf orbitals by neutrons. Hence electron 
captures on such nuclei are only possible due to correlations of nucleons across the 77 = 40 shell gap. 
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Such cross-shell correlations are a rather slowly converging process and usually require multi-particle- 
multi-hole excitations [31]. Taking this into account, the electron capture rates for heavier nuclei have 
been calculated within a hybrid approach. At hrst, the Shell Monte Carlo Approach (SMMC) [32, 33] 
has been used to determine occupation numbers of proton and neutron orbitals by performing SMMC 
calculations for individual nuclei at hnite temperature within the complete pf-sdg shell and taking 
the appropriate nuclear correlations into account. Using these hnite-temperature occupation numbers 
electron captures have been calculated using the RPA approach. The latter step is justihed as at 
the conditions under which the nuclei with A > 65 are relevant the electron chemical potential is 
noticeably larger than the respective nuclear Q values making the capture rates mainly sensitive to the 
centroid and total value of the transition strength, but not to its detailed structure [28]. Forbidden 
transitions contribute increasingly to the capture rates with growing density and, relatedly, electron 
chemical potential. Inspired by the IPM supernova simulations have for a long time assumed that 
electron captures on heavy nuclei (with N > 38 [34]) vanish and that the capture occurs solely on free 
protons (which have the larger individual capture rate than heavy neutron-rich nuclei due to the smaller 
Q value, but are much less abundant in the core caused by the low entropy). More recent supernova 
simulations have incorporated the modern shell model rates and hnd that the capture on heavy nuclei 
dominates during the entire collapse [35, 36, 37]. By measuring the GT_|_ strength for ^®Se (with 34 
protons and 42 neutron) it has been experimentally demonstrated that the strength is indeed unblocked 
by nuclear correlations across the A^ = 40 gap [38]. Data from transfer reaction measurements [39], in 
agreement with large-scale shell model calculations [40], indicate about 4 neutron holes in the pf shell 
in the ^®Se ground state, unblocking GT+ transitions of valence protons. 

Electron capture on nuclei is the main neutrino source during collapse. Fig. 1 shows the respective 
neutrino spectra produced at different moments during the core collapse of a 15 Mq and a 25 Mq star, 
respectively. The spectra have been calculated from the individual electron captures of about 3000 
nuclei, averaged over an appropriate NSE distribution. The collapse snapshots are chosen to support 
the following discussion on neutrino trapping and thermalization. Due to the hnite temperature of the 
environment and the presence of many nuclei the neutrino spectra are quite wide. At densities of a few 
10^^ g cm“^ the spectrum is centered around neutrino energies of order 10 MeV. Due to the increase 
of electron chemical potential with density (~ the neutrino spectra are shifted to higher energy 
with growing density. At p > 10^^ g cm“^ the centroid is at neutrino energies E^, cs 20 MeV. Screening 
due to the astrophysical environment effectively reduces the chemical electron potential and enhances 
the nuclear Q value. These effects decrease the electron capture rates slightly and shift the neutrino 
spectra to lower energies (by about 2 MeV as is shown in Fig. 1). 

An important change in the physics of the collapse occurs, when neutrinos start to get trapped in 
the core at densities ptrap ~ 4 x 10^^ g cm“^. The main process is elastic scattering of neutrinos on 
nuclei. The respective total cross section for this process is given by [3] 

aei(E,) = ^ [N-{l-4sm^9w)zYEl (1) 

where E^ is the energy of the incident neutrino; Gp and 6w are the Fermi coupling constant and the 
Weinberg angle (sin^ ~ 0.231), respectively. As (1 — 4sin^6*w) ~ 0.08, elastic neutrino scattering 
occurs mainly as a coherent process on the neutrons in the nuclei. In Eq. 1 it has been assumed that 
the elastic form factors for protons and neutrons can be replaced by Ez,N{,<f) = Ez^n{q‘^ = 0) = 1 
which is justihed at the momentum transfers involved. For odd-A and odd-odd nuclei with ground 
state angular momenta J ^ 0 and at hnite temperature there can also be a GT contribution to the 
cross section. However, due to the strong fragmentation of the GT strength this is usually negligibly 
small. Furthermore under supernova conditions even-even nuclei dominate the nuclear composition. As 
is demonstrated in [3] neutrinos with energies about 20 MeV will have a mean free path of order 0.5 km 
at densities 10^^ g cm“^ which is signihcantly smaller than the core radius at the same density. As a 
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Figure 1: Neutrino spectra produced by electron captures on nuclei at two different densities (3.76 x 
10^^ g cm“^ and 2.47 x 10^^ g cm“^) during the collapse of a star with 15 Mq and at 1.13 x 10^^ g cm“^ 
and 1.54 x 10^^ g cm“^ for a 25 Mq star. For the solid line an NSE distribution has been assumed 
for the matter composition; the dashed line shows the spectrum if heavy nuclei (with A > 100) are 
neglected in the distribution. The dotted curve shows the spectra if screening corrections due to the 
astrophysical environment are not included, (from [28]) 
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Figure 2: Comparison of the central values for the electron (Yg) and lepton Y^ep abundances (left panel), 
temperature and entropy (middle panel) and electron neutrino emission rate (right panel) obtained in 
a spherical simulation of the collapse of an 18 Mq star using the shell model electron capture rates 
(LMSH) or neglecting capture on heavy nuclei with N > 38 (Bruenn rates [34]). The insert in the right 
panel compares the average neutrino energies, (courtesy Hans-Thomas Janka, adapted from [45]) 


consequence neutrinos scatter quite often in the core and their propagation has to be described by a 
diffusion process. Indeed, the diffusion time scale at p = 10^^ g cm“^ is longer than the competing 
collapse time scale and neutrinos are effectively trapped during the hnal phase of the collapse. 

Neutrinos exchange energy with the core matter by inelastic scattering, mainly on electrons. As 
electrons are highly degenerate at this stage of the collapse, they can only gain energy and hence 
neutrinos are down-scattered in energy. By inelastic scattering on nuclei neutrinos can both loose 
and gain energy (as is discussed in Section 4). Simulations, however, show that this process is less 
important during collapse than inelastic neutrino-electron scattering which is found to be quite efficient 
in thermalizing the neutrinos with the rest of the matter. This is effectively achieved at densities in 
excess of 10^^ g cm“^ as is demonstrated in Fig. 2 which shows results from a (spherical) collapse 
simulation of an 18 Mq star. The hgure shows that the lepton abundance, i.e. the sum of abundances 
of electrons and electron neutrinos, reach a constant value once neutrinos are trapped and a Fermi sea of 
electron neutrinos builds up in the core. The hgure also clearly demonstrates the impact of the modern 
electron capture rates (LMSH) on the collapse dynamics if compared to results based mainly on the 
IPM approach (neglecting electron capture on nuclei with neutron numbers N > 38 [34]). As electron 
captures on nuclei dominate over those on free protons the electron neutrino emission rate is larger, 
while the central core values for temperature, entropy and are smaller. At densities p < 3 x 10^° 
g cm“^ nuclei from the mass region A = 45 — 65 dominate. For these nuclei the electron capture rates 
obtained by shell model diagonalization [27] are noticeably smaller than those inspired by the IPM (and 
corrected by data whenever available) [41, 42, 43, 44]. As neutron-rich nuclei have larger Q values than 
protons, the average neutrino energies are smaller in simulations incorporating the shell model capture 
rates. 

After neutrino trapping, the collapse proceeds homologously [46], until nuclear densities = 













2.5 X 10^"^ g cm“^) are reached. Since nuclear matter has a much lower compressibility, the homologous 
core decelerates and bounces in response to the increased nuclear matter pressure. This drives a shock 
wave into the outer core; i.e. the region of the iron core which lies outside the homologous core and in the 
meantime has continued to fall inwards at supersonic speed. The core bounce with the formation of a 
shock wave is the mechanism that ultimatively triggers a supernova explosion, but the exact mechanism 
of this physically appealing scenario are still uncertain and controversial. It appears as if the energy 
available to the shock is not sufficient, and the shock uses up its energy in the outer core mostly by 
the dissociation of nuclei into nucleons [3, 37]. This change in composition results to additional energy 
losses, because the electron capture rate on free protons is signihcantly larger than on neutron-rich 
nuclei due to the higher Q-values of the latter. A large fraction of the neutrinos produced by these 
capture behind the shock leave the star quickly in what is called the neutrino burst at shock break out, 
carrying away energy. This leads to further neutronization of the matter. The shock wave is weakened 
so much that it hnally stalls and turns into an accretion shock. 

After the core bounce, a compact remnant begins to form in the center. Depending on the stellar 
mass, this is either a neutron star (for progenitor stars with masses roughly smaller than 25 Mq) or 
a black hole. The nascent proto-neutron star contains a large number of degenerate electrons and 
neutrinos, the latter being trapped as their mean free paths in the dense matter is significantly shorter 
than the radius of the neutron star. As the trapped neutrinos diffuse out, they convert most of their 
initially high degeneracy energy to thermal energy of the stellar medium [47]. The cooling of the 
protoneutron star then proceeds by pair production of neutrinos of all three generations which diffuse 
out. After several tens of seconds the star becomes transparent to neutrinos and the neutrino luminosity 
drops signihcantly [48]. 

In what is called the ‘delayed neutrino-heating mechanism’ the stalled shock wave can be revived 
by neutrinos [49]. These carry most of the energy set free in the gravitational collapse of the core [47] 
and deposit some of it in the layers between the nascent neutron star and the stalled prompt shock, 
mainly by absorption on nucleons. This lasts for a few 100 ms, and requires about 1% of the neutrino 
energy to be converted into nuclear kinetic energy. The energy deposition increases the pressure behind 
the shock. However, it appears that energy transport by neutrinos is not enough to explode models 
with spherical symmetry [50, 51, 52]. However, multi-dimensional hydrodynamical simulations in two 
dimensions [53, 54, 55, 56] and in 3 dimensions [57, 58, 59, 60, 61] show that the region of low density, but 
rather high temperature (named ’hot neutrino bubble’), which forms between the neutron star surface 
and the shock front, is convectively unstable developing large convective overturns which enhance the 
neutrino energy deposit behind the shock [54, 62, 63, 64]. This is additionally aided by instabilities 
of the shock to non-radial deformation (standing accretion shock instability SASI) which increases the 
time matter stays in the heating region [65, 66, 67, 68, 69, 70]. 

The simulations sensitively depend on the nuclear Equation of State (EoS) which is still not suffi¬ 
ciently well-known. It strongly influences the formation, the strength and the evolution of the shock 
wave, the latter indirectly via the compactness and the contraction of the nascent neutron star. For 
densities less than about pat/IO, nucleonic matter exists of individual nuclei. At larger densities matter 
exists in a variety of complex shapes caused by a competition of repulsive Coulomb forces and nuclear 
attraction until it reaches uniform nuclear matter. There are several EoS [71, 72, 73, 74, 75, 76, 77] 
frequently used in supernova simulations. A novel ansatz, based on the relativistic mean-field with 
density-dependent couplings and the explicit consideration of light-particle bound states, has recently 
been presented [78]. Refs. [79, 80] has studied the implications of an Equation of State with an explicit 
QCD phase transition to quark matter on the postbounce supernova evolution. 

In the dense environment encountered in the late stage of the collapse or in the nascent neutron 
star, neutrino processes are modified by ion-ion correlations [81]. This includes the coherent elastic 
scattering of neutrinos on nuclei [82, 83, 84]. Correlations are particularly important for neutrino-pair 
bremsstrahlung, NN ^ NNvv. Calculations performed on the basis of Chiral Effective Field Theory 
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Figure 3: Results from 2D core-collapse simulations with two different equations of states (“Wolff” [73] 
and “L&S” [71]) for a 15 Mq star, using the neutrino transport technique of Refs. [89, 90]. The figure 
displays the neutrino luminosities with the prompt Vf. burst in the left window and the postbounce 
luminosities of z/g (solid lines), he (dotted) and heavy-lepton z/’s and z/’s (dashed) in the right window. 
The results of the 2D runs were averaged over latitudinal angles 9. (from [37]) 


have recently shown that noncentral contributions to the strong interaction, i.e. tensor and spin-orbit 
forces, reduce the neutrino pair production rate significantly compared to one-pion exchange approx¬ 
imations usually adopted in supernova simulations [85, 86]. Neutrino-pair bremsstrahlung, together 
with neutrino-pair annihilation into pairs of other flavor (mainly z/ghe —?• Vx^xi x = is crucial 

for transport of fi and r neutrinos [87, 88] and its inclusion in supernova simulations has resulted in 
noticeable changes of the spectrum of these two neutrino types towards smaller average energies. 

Fig. 3 shows the time evolution of the luminosities in the different neutrino flavors obtained in a two- 
dimensional simulation for two different Equation of States. During the collapse phase the supernova 
overwhelmingly emits electron neutrinos, produced by electron capture on nuclei (the production rate 
of other neutrino types by nuclear deexcitation during collapse is much less and is discussed in Section 
4). The sharp reduction in the z/g luminosity just before bounce (at tpb = 0) is related to neutrino 
trapping in the last phase of the collapse. It is followed by the sharp peak of the neutrino burst just 
after bounce, caused by photodissociation of nuclei into free nucleons and the associated fast electron 
captures on free protons. After about 50 ms the neutrino luminosities are mainly due to neutrino pair 
production caused by electron-positron annihilation in the dense and hot environment. Consequently 
supernova shine now in all neutrino flavors. As /i and r neutrino pairs can couple to pairs only 
by neutral current, their production rate is smaller than the one of z/gZ/g pairs. 

The radius from which neutrinos can be considered to be Tree streaming’ is called the neutri- 
nosphere [3]. For electron neutrinos the neutrinosphere approximately corresponds to the region with 
trapping densities of a few 10^^ g cm“^, which is at radii of about 20 km in the later postbounce 
phase [91]. As neutrino cross sections depend on the neutrino energy, the neutrinosphere is ’diffuse’; 
i.e. its radius decreases with neutrino energy. 
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On quite general grounds one expects a hierarchy in the average energies of neutrinos of different 
flavors: As fi and r neutrinos have not sufficient energies to interact with matter by charged-current 
reactions, they decouple at smaller radii, and higher temperatures, than electron neutrinos. Further¬ 
more, as the matter is neutron rich, anti-electron neutrinos, which by charged-current only interact with 
protons, decouple at smaller radii and higher temperature than electron neutrinos, which interact by 
charged-current reactions with neutrons. The hierarchy in temperatures at which the neutrino flavors 
decouple translates then into a hierarchy of average energies: > {EpJ > {E^J. For current 

estimates for the average neutrino energies see refs. [6, 8]. 

It is one of the major aims for future supernova neutrino detections to conhrm this hierarchy. 
However, these measurements will depend strongly on neutrino flavor oscillations along the way to the 
detector. One possibility are adiabatic flavor conversions [92] which implies noticeable mixing of the Ug 
and z/a, spectra in case of the normal mass hierarchy and a complete swap for the inverted mass hierarchy 
[92]. For our nucleosynthesis discussion in Section 5 the flavor conversion induced by neutrino-neutrino 
refraction [93] is more important. In fact, several studies [94, 95, 96, 97] indicate that collective neutrino 
flavor oscillations, which occur in the high-neutrino-density environment surrounding the neutron star, 
swap the spectra of z/g and neutrinos in certain energy intervals bounded by sharp spectral splits. 
This can have interesting consequences for the nucleosynthesis in the up process (see Section 5), 

3 Theoretical background and experimental validation 

In this Section we want to discuss the theoretical background to describe neutrino-induced reactions on 
nuclei at neutrino energies relevant for supernova dynamics and nucleosynthesis. We will at hrst focus 
on charged- and neutral-current reactions on nuclei in their ground state. These reactions are directly 
relevant for studies of the neutrino-induced response of detectors for supernova neutrinos. The hnite 
temperature in the supernova environment implies that nuclei exist in an ensemble with excited nuclear 
states thermally populated. This requires special attention. As charged-current (z/e,e“) and (z/g, e"*") 
reactions are the inverse of stellar electron and positron captures, the relevant reaction rates can be 
obtained from the latter by applying the principle of detailed balance. In fact, in this way (z/g, e~) 
and (z/g, e"*") reactions on nuclei are included in supernova simulations for a long time. Furthermore, 
charged-current reactions induced by muon and tan neutrinos are neglected in supernova simulations 
as their energies are not sufficient to initiate these processes at supernova conditions where individual 
nuclei exist. In the supernova neutral-current reactions on nuclei might alter the neutrino spectrum 
and/or contribute to the energy exchange between matter and neutrinos. Hence hnite temperature 
effects are relevant and we will discuss below how these can be estimated and which inhuence they have 
on cross sections. The neutrino energies are often large enough to excite the nucleus to states above 
particle thresholds. In the hot supernova environment the subsequent particle decay can be neglected as 
the matter composition can be described in NSE. However, the decay is important for nucleosynthesis 
studies at lower temperatures as it alters the abundance distributions. The relevant decay probabilities 
will be described within the framework of the statistical model. 

This Section is divided in 3 subsections in which we will derive the relevant cross sections, summarize 
briefly the hierarchy of nuclear models adopted to describe the nuclear structure problem at the various 
energy scales involved in the neutrino-nucleus reactions of interest and hnally validate this ansatz by 
comparison to relevant experimental data. 

3.1 Cross sections 

For the semileptonic processes 

+ Z^N —t z+l^*N-l + e ; t'g + zAjv —t z+1^*N-1 + S"*") ^ + Z^N —t z^*N + (2) 


11 


the excitation cross section from the ground states of a nucleus (dehned by proton and neutron numbers 
Z and respectively) to a discrete hnal state is given in the extreme relativistic limit (hnal lepton 
energy e/ » lepton mass rrifC^) by [98, 99]: 



(dai^f\ Gl-e) dvrcos^f 

oo oo 

<^T 

.,7=0 J=1 


(3) 

where 

ai,p = \{Jf\\Mj{q) + -^LMU)\^ 


(4) 

and 

^T= (-§^+lan2f) X [|(J/||J7“^(g)||Ji)P + |(J/||Jj(g)||Ji)P] 

T2tanf y/^ + tan7 X Re [(J/Pr‘’(9)ini>V/ll7(9)in<>1 ■ 

(5) 


Here Gp is the Fermi coupling constant, 0 the angle between the incoming and outgoing lepton, and 
= (a;, ^ (<? = |^) the four-momentum transfer. The minus-(plus) sign in Eq. (4) refers to the neutrino 
(antineutrino) cross section. The quantities Mj, L j, Jf and denote the multipole operators for 
the charge, and the longitudinal and the transverse electric and magnetic parts of the four-current, 
respectively. Following Ref. [98] they can be written in terms of one-body operators in the nuclear 
many-body Hilbert space. The cross section involves the reduced matrix elements of these operators 
between the initial state Jj and the hnal state Jf. 

The Fermi function F{Z, e/) appears in the cross section for the charged-current process and accounts 
for the Coulomb interaction between the hnal charged lepton and the residual nucleus. It can be derived 
relativistically as the Coulomb correction obtained from a numerical solution of the Dirac equation for 
an extended nuclear charge [ 100 , 101 ]: 


e/) • Lq 

F(27 + 1) 

where Z denotes the atomic number of the residual nucleus in the hnal channel, ej the total lepton 
energy (in units of rrtic^) and pi the lepton momentum (in units of m;c), R is the nuclear radius (in 
units of A = h/{mic)) and 7 and y are given by (a = hne structure constant): 

7 = y/l — (aZ)"^ and y = aZ — . (7) 

The numerical factor Lq in (5), which describes the hnite charge distribution and screening corrections, 
is nearly constant(~ 1 . 0 ), and can be well approximated by a weakly decreasing linear function in 
Pi [101]. 

At higher energies, the Fermi function for s waves becomes a poor approximation when higher partial 
waves also contribute for piR > 1 . Guided by the distorted-wave approximation of quasielastic electron 
scattering, in that case the Coulomb ehects can be treated in the ‘Ehective Momentum Approximation’ 
in which the outgoing lepton momentum pi is replaced by the ehective momentum 


with Fo{Z, ej) = 4{2piR)‘^^^-^^ 


PeS = \l -E'eff “ = Ei - Vc{R) , (8) 

where Vc{R) = 3e^Z/(2i?) is the Coulomb potential at the origin. In the formalism above [102], the 
Coulomb ehect is taken into account not only by using the ehective momentum, but also by replacing 
the phase space factor piEi by Peff-Ees (see also [103] where this procedure is called modihed ehective 
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momentum approximation, and shown to work quite well). In practice, a smooth interpolation between 
these two regimes of treatment of the Coulomb effects is used. 

The total cross section is obtained from the differential cross sections by summing over all possible 
hnal nuclear states and by numerical integration over the solid angle dVL. 

Supernova neutrino energies are sufficiently high to excite nuclei in the hnal channel to energies above 
particle thresholds. The nuclear state will then subsequently decay by particle emission. Under certain 
conditions the hnal nucleus might even be excited to energies high enough to allow for multi-particle 
emissions. Such a situation occurs in neutral-current reactions induced by supernova neutrinos as 
their expected average energies of about 18 MeV often suffices to excite nuclei above the 2-particle 
emission threshold. Such a decay is also relevant for charged-current reactions induced by supernova 
Vf. neutrinos if the Q value is large enough. Such a situation will certainly occur for r-process nuclei if 
they are subject to charged-current reactions with supernova Ve neutrinos. As nuclei on the r-process 
path also have very small neutron separation energies, the state excited by the neutrino process in the 
daughter nucleus will decay by the emission of several neutrons [104, 105]. It has been pointed out that 
neutrino-induced hssion of r-process nuclei will also be accompanied by multiple neutrons [106, 107]. 

Neutrino-induced reactions with subsequent particle emission can be described in a two-step process. 
In the hrst step a nuclear level / at energy Ef with angular momentum Jf and parity vtj is excited. The 
respective reaction cross section is described by Eq. 3. The subsequent particle decay is calculated on the 
basis of the statistical model. At relatively low excitation energies, where one-particle decays dominate, 
these calculations are preferably done on the basis of statistical model codes which explicitly consider 
experimentally known nuclear spectra. There are several such codes existing: SMOKER [108], MOD- 
SMOKER and SMARAGD [109], TALYS [110]. There have been a few applications where such codes 
have been extended to treat also the cascade of two-particle emissions. However, these treatments are 
computationally quite demanding, thus multi-particle emissions are usually studied within a Monte- 
Carlo treatment of the decay cascades. Such a dynamical code is ABLA07 [111] that describes de¬ 
excitation of the compound nucleus through the evaporation of light particles (protons, neutrons, a and 
7 ) and hssion. The particle evaporation is considered in the framework of the Weisskopf formalism [112], 
while in the description of hssion the ABLA code treats explicitly the relaxation process in deformation 
space and the resulting time-dependent hssion width using an analytical approximation [113] to the 
solution of the Fokker-Planck equation. The wide range of excitation energies and the large variety 
of nuclei of interest demands for a consistent treatment of level densities as a function of excitation 
energy and deformation. The nuclear level density is calculated according to the Fermi-gas formula 
taking into account the shell corrections, collective excitations and pairing correlations; details on the 
implementation of these ehects in the nuclear level density, their damping with excitation energy and 
their deformation dependence are described in details in Refs. [114, 115, 116]. The angular-momentum- 
dependent hssion barriers are taken from the hnite-range liquid-drop model [117]. In order to describe 
the hssion-fragment mass and charge distributions, as well as their kinetic energies, the ABLA code is 
coupled to the semi-empirical Monte-Carlo hssion model PROFI [118]. 

3.2 Hybrid model for allowed and forbidden contributions 

Reactions of low-energy neutrinos with < 20 MeV oh nuclei are dominated by allowed J = O’*" and 
J = transitions (see Fig. 4). In the limit of vanishing momentum transfers q = 0, which is justihed 
at the low neutrino energies involved, these two transitions between initial {i) and hnal (/) nuclear 
states reduce to the Fermi and Camow-Teller transitions, respectively: 


Bt,(F) = T(T+l)-T,,T,,, (9) 

where the transition is only possible between Isobaric Analog states and where we have neglected the 
reduction in the overlap between nuclear wave functions due to isospin mixing which is estimated to be 
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Figure 4: Multipole decomposition of (z/g, e~) cross 
(from [119]) 



sections on ®®Fe for two different neutrino energies. 


small 0.5% [120]), and 


B,f(GT) 


\9v / eff 2 Jj + 1 


( 10 ) 


where the matrix element is reduced with respect to the spin operator cr only (Racah convention [121]) 
and the sum runs over all nucleons. For the isospin rising and lowering operators, t± = :^iTy)/2, we 

use the convention t+p = n; thus, ‘+’ refers to electron capture and f3~^ transitions and ‘ — ’to positron 
capture and (3~ transitions. Finally, {gA/gv)eS is the effective ratio of axial and vector coupling constants 
that takes into account the observed quenching of the GT strength [122]. One often uses [123, 124, 125] 


9v 


— ] = 0.74 I — 


eff 


9V 


( 11 ) 


bare 


with {qa/ 9v)ha.re = —1.2599(25) [120]. If the parent nucleus (with isospin T) has a neutron excess, then 
the GT_ operator can connect to states with isospin T —1, T, T + 1 in the daughter, while GT_|_ can only 
reach states with T + 1. This isospin selection is one reason why the GT_|_ strength is more concentrated 
in the daughter nucleus (usually within a few MeV around the centroid of the GT resonance), while 
the GT_ is spread over 10-15 MeV in the daughter nucleus and is signihcantly more structured. 

The GT response is in turn sensitive to nuclear structure effects and hence its contribution to the 
cross section has to be derived from a model which is capable to describe the relevant nuclear structure 
and correlation effects. This model is the diagonalization shell model [31]. Higher multipoles will, 
however, contribute to the cross section for larger neutrino energies. For each of these multipoles the 
response of the operator will be fragmented over many nuclear states. However, most of the strength 
resides in a collective excitation, the giant resonance, whose centroid energy grows with increasing 
rank A roughly like Xhuj ~ 41A/H^/^ MeV. As the phase space prefers larger hnal neutrino energies, 
the average nuclear excitation energy grows noticeably slower than the initial neutrino energy. As a 
consequence, initial neutrino energies are noticeably larger than the energy of a giant resonance, when 
the later will contribute to the neutrino-nucleus cross section. Fortunately, the neutrino-nucleus cross 
section depends then mainly on the total strength of the multipole excitation and its centroid energy, 
and not on the detailed energy distribution of the strength (as this is the case for the Gamow-Teller 
response at low neutrino energies). Thus the higher multipole contributions to the neutrino-nucleus 
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cross section can be derived within the RPA making nse of the fact that the RPA describes the energy 
centroid and total strength of mnltipoles other than the Gamow-Teller qnite well. 

In snch a hybrid model, the total reaction cross section consists then of two parts 

= + ( 12 ) 

where the shell model part, cr^™, describes the allowed contribntions to the cross section for (z/g, e“) 
(GT_) and (Pe,e+) (GT+) reactions and is given by: 

^ k{ElF(±Z ± 1, El) (,B,(F) + B/(GT±)) (13) 

f 

where Gp is the Fermi constant, 6c the Gabibbo angle, and and Ef, the momentnm and energy of the 
ontgoing electron or positron. The snm is over the final (/) nnclear states with energies Ej. There is no 
Fermi contribntion to (z/g, e"*") reactions on nuclei with neutron excess. For inelastic neutrino scattering 
to an excited nuclear state at energy Ef the allowed contribution to the cross section is given by: 

^ El, jB,(G%) (14) 

TT ^^ 

/ 

where energy conservation yields E^j = E^ — Ef. These allowed contributions replace the A = O’*" and 
A = !■*■ contributions in the cross section formula above. All other multipole contributions can be derived 
adopting the RPA as the nuclear model to describe the initial and hnal nuclear states. We mention that, 
for light nuclei like ^^G, shell model calculations are now also possible for the first-forbidden transitions 
(A = 0“, 1“, 2“) requiring model spaces spanned by 3hu [126, 127, 128]. 

3.3 Inelastic neutrino-nucleus reactions at finite temperature 

So far, our formalism assumes that the target nucleus is in the ground state corresponding to an 
environment with temperature T = 0. This is, of course, appropriate for neutrino-induced reactions 
in earthbound detectors and it is a reasonable approximation for the supernova environment in which 
neutrino-nucleus reactions occur in nucleosynthesis processes (T < 100 keV). However, the T = 0 
approximation is not valid if we want to study inelastic neutrino-nucleus reactions in the hot supernova 
environment (with temperatures T > 1 MeV). (We recall that charged-current reactions under these 
conditions can be derived from electron and positron captures by detailed balance.) Here nuclei exist as 
a thermal ensemble with excited nuclear states at energy Ei with spin Jj populated with a Boltzmann 
weight of Gi = (2Jj -|- 1) exp{—Ei/kT). There are two interesting effects which occur at finite T: i) 
neutrinos can be up-scattered in energy; i.e. the neutrino picks up energy from the nucleus which 
by the scattering reaction deexcites; ii) at zero temperature for even-even nuclei with their J = O’*" 
ground states GT transitions can only be induced by neutrinos with minimum energies to bridge the 
gap to the first excited J = 1^ state in the nucleus. Transitions to excited states at lower excitation 
energies can be induced by higher multipoles, but the respective cross sections are usually quite small. 
As a consequence inelastic neutrino-nucleus cross sections for even-even nuclei drop quite fast with 
decreasing small neutrino energies. This is not expected to be the case for inelastic neutrino-nucleus 
reactions at hnite temperatures. Due to the exponential increase in the nuclear level density states 
at modest excitation energies are likely to have states in their energy vicinity with which they can be 
connected by GT transitions. Furthermore, it is even likely that GT transitions (or with less strength 
by other multipoles) can lead to levels at lower energies, a transition which is in principle possible for 
all neutrino energies. In summary, we expect on general grounds that the cross section is in particular 
at low neutrino energies modified at finite temperature. 
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Unfortunately an evaluation of the inelastic cross section at finite temperature by explicitly summing 
over the thermally populated states is computationally unfeasible due to the many states involved at 
temperatures of order MeV. There are two suggestions how to overcome this problem and treat inelastic 
neutrino scattering at hnite temperature. We briefly discuss both in turn. 

Refs. [129, 130] apply Brink’s hypothesis which states that for a given excited nuclear level i the 
strength distribution built on this state, Si{E), is the same as for the ground state, So{E), but shifted 
by the excitation energy Ef. Si{E) = Sq{E — Ei). The validation of Brink’s hypothesis has been in 
detailed studied in ref. [131], however, for electron capture on sd-shell nuclei (see also [27]). Upon 
applying the Brink hypothesis, the down-scattering part becomes independent of temperature and can 
be solely derived from the ground state GT distribution. With this approximation, the Gamow-Teller 
(shell model) contribution to the cross section becomes: 


cr.. 
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TT 
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■ f f 


(15) 


where Gp is the Fermi constant, and E^/^if is the energy of the scattered neutrino, Ey/^if = Ey + {Ei—Ef), 
with Ei, Ef denoting the initial and hnal nuclear energies. 

The hrst term in Eq. (15) arises from Brink’s hypothesis. By construction, this Brink term does not 
allow for neutrino up-scattering. These contributions to the cross section are comprised in the second 
term, where the sum runs over both initial {i) and final states (/). The former have a thermal weight of 
Gi = (2Jj -|- 1) exp{—Ei/kT), where Ji is the angular momentum, Ei is energy of the initial state, and 
T is temperature; G = Gi is the nuclear partition function. The up-scattering contributions are the 
more important i) the lower the nuclear excitation energy (Boltzmann weight in the thermal ensemble), 
ii) the larger the GT transition strength Bif, and iii) the larger the final neutrino energy. Guided by 
these general considerations Ref. [130] has approximated the second term in Eq. (15) by explicitly 
considering GT transitions between nuclear states with Ei > Ef, where the sum over the hnal nuclear 
states is restricted to the lowest nuclear levels. As is usual in astrophysical applications, experimental 
data (in this case mainly excitation energies) have been used whenever available. 

As stressed above, the derivation of the hrst (down-scattering) term in the cross section explicitly 
assumes that the distributions on excited states is the same as for the ground state. This is likely not 
the case as the vanishing of pairing correlations with increasing temperature should ehect the relative 
excitation energy of the strength centroid and move it slightly down in energy. In fact Shell Model 
Monte Garlo studies, performed at hnite temperature, conhrm this conjecture [132]. To account for 
this ehect, Dzhioev et al. developed a model which extends the QRPA approach to hnite temperature 
(Thermal Quasiparticle Random-Phase Approximation, TQRPA [133, 134]). In the applications the 
authors consider multipole contributions up to A = 3. However, it is illustrative to assume that the 
scattering process is given solely by the contribution from the allowed J = 1+ multipole. Then the 
cross section induced by neutrinos with energy Ey at temperature T can be written as 

a{Ey, T) = ^ V(E, - Uif^i + ^ y^{Ey + uf exp(-^)<F,. (16) 

TT TT Kl 


Analog to Eq. (15) the hrst term corresponds to the down-scattering, and the second to the up- 
scattering contributions. The sum in the down-scattering term runs over all I"*" thermal phonon states 
with positive energy Ui < Ey. The Boltzmann factor suppresses the contribution from highly thermally 
populated states in the up-scattering term. These transitions are, however, somewhat favored by the 
phase factor (Ey + u)'^ [134]. In Eq. (16) denotes the transition strength of the GT operator. We note 
that the TQRPA approach obeys detailed balance which, however, is slightly violated in the approach 
based on Eq. (15). 

Using iron isotopes as example we will discuss results from both approaches in Section 4. 
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3.4 Nuclear models and their validation by data 


Measurements of neutrino-nucleus cross sections are very rare. For charged-current reactions there 
exist the set of measurements for solar electron neutrinos on the nuclei comprising the various solar 
neutrino detectors. These data - despite being very valuable for solar physics and the proof of neutrino 
oscillations - are only of restricted use to check the applicability of nuclear models to describe neutrino- 
nucleus reactions for supernova neutrinos. Solar neutrinos have an energy spectrum, which is well 
understood from the reactions mediated by the weak interaction in solar hydrogen burning and distorted 
by matter-enhanced oscillations in the Sun. However, solar neutrino energies are too low to excite giant 
nuclear resonances which in turn are expected to dominate the cross sections for supernova neutrinos. 
This means that transitions induced by solar neutrinos are strongly state-dependent and are quite a 
challenge to model. There have been a few theoretical studies performed for solar neutrinos within the 
framework of effective held theory [135], shell model and RPA. 

For the deuteron, which has been the target nucleus in the SNO detector, neutrino-scattering calcu¬ 
lations have been performed in effective held theory (EFT). This approach allows to split the neutrino- 
induced cross section into two pieces cr{E) = a{E) + where a{E) and b{E) are known and 

all unknown ehects can be lumped together in one unknown parameter Lia (isovector two-body axial 
current) which must be determined by experiment [136, 137]. Although a{E) contributes dominantly to 
the cross section, the current uncertainty in Lia (about 50%) is still large [138, 139, 140, 141]. Similar 
results have been obtained based on the phenomenological Lagrangian approach [142] 

For the two nuclei and ®®Fe there exist measurements of cross sections induced by z/g neutrinos 
[143, 144, 145, 146, 147, 148]. The neutrinos are produced by muons decaying at rest and have therefore 
spectra which are slightly shifted to larger neutrino energies than supernova neutrinos. Nevertheless the 
respective neutrino-nucleus cross sections are expected to be dominated by contributions from allowed 
(Fermi and GT) and hrst-forbidden transitions and are therefore valuable test cases for nuclear models. 
In turn measurements of these transition strengths serve as strong indirect constraints of the ability of 
nuclear models to describe neutrino-nucleus scattering at supernova neutrino energies. 

In Ref. [18] it is argued on the basis of general considerations that, depending on the neutrino 
energies, the nuclear structure problem has to be treated with different demands and details. As we 
will show below these requirements are fulhlled by a hierarchy of models which is actually the basis 
of the hybrid model dehned above: the nuclear shell model for low neutrino energies and the Random 
Phase Approximation (or its variants) for higher energies. 

i) Low energy neutrino scattering requires the detailed reproduction of the nuclear response (mainly 
allowed transitions) at low nuclear excitation energy. The method of choice here is the shell model which, 
for a hxed model space of valence nucleons, accounts for nucleon-nucleon correlations among these 
valence particles via an effective interaction. Shell model wave functions have good quantum numbers 
for angular momentum, parity and isospin. Modern shell model codes allow now diagonalization in 
model spaces with up to a few billion conhgurations. Based on a Monte Carlo strategy of importance 
sampling shell model calculations for low energy spectra have been extended to even larger model spaces. 
Hence shell model diagonalization is now possible for medium-mass nuclei in the iron-nickel region in the 
full pf shell as is desirable for calculations of GT transitions.. For lighter nuclei, like and shell 
model calculations can be performed in multi-ho; model spaces allowing also for fully correlated studies 
of hrst-forbidden transitions. Such transition strengths are calculated employing the Lanczos technique 
for a hnite number (of order 100) iterations. As a consequence the lowest transitions correspond to well- 
dehned (shell model) states, while at higher excitation energies the calculated strength corresponds to 
the summed strength per energy interval (which depends on the number of iterations considered). For 
a review on the shell model and its abilities and applications the reader is refered to Ref. [31]. 

ii) The RPA has been developed to describe the collective excitations of a nucleus by considering 
the one-particle one-hole excitations of the correlated ground state. Gompared to the shell model. 
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the RPA calculations usually use model spaces with a (signihcantly) larger set of single particle, but 
they consider, however, only a strongly truncated part of the possible conhgurations. This has the 
consequence that the strength is less fragmented in RPA than in shell model calculations. 

It is important that in the shell model and in the RPA the total strength is hxed by sum rules. Of 
particular relevance for the description of neutrino-nucleus reactions at supernova neutrino energies are 
the Ikeda sum rule for the Gamow-Teller strength 


GT{Z ^ Z + GT{Z ^ Z-l)i = 3{N - Z), 


and the Thomas-Reiche-Kuhn sum rule for the dipole strength 


— Eq)B(E1; 0 —>■ f) 
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where Mn is the nucleon mass. For nuclei with large neutron excess, the Ikeda sum rule basically hxes 
the total GT_ strength relevant for reactions. 

It is known from studies of (3 decays and of GT distributions [123, 125] that both shell model 
and RPA calculations systematically overestimate the GT transition strengths. For the shell model 
calculations this shortcoming can be overcome by introducing a constant quenching factor via the 
introduction of an effective axial-vector coupling constant which for pf shell nuclei has empirically been 
determined as g'^ = O.TAgA [125]. GT strengths, obtained in RPA calculations, are often quenched by 
the same factor, although this procedure in RPA studies is hampered by the limited consideration of 
conhgurations within the model space. 

Replacing the full A = !■*■ operator by its approximation in the q = 0 limit (the GT operator) is 
only justihed for small incident neutrino energies. At modest neutrino energies the consideration of 
the hnite momentum transfer reduces the cross section. In Ref. [149] it is suggested to correct cross 
sections calculated from the shell model GT distribution by the ratio of RPA calculations performed 
for hnite momentum transfer and for g = 0, respectively. For neutrinos produced by muon decay at 
rest the correction is of order 20%. For supernova neutrino energies the correction is smaller. 

Ref. [18] gives a detailed comparison of data obtained for charged-current reactions on ^^G induced 
by muon-decay-at-rest (DAR) z/g neutrinos with various model calculations. The exclusive reaction to 
the J = 1^ ground state of is given by a GT transition. In general models agree quite well with 
experiment (see Table I in Ref. [18]). However, it must be noted that the calculations are usually 
constrained by other experimental data which hx the GT matrix element (like the (3 decays of the 
and ground states or the muon capture to the ground state). Hence these calculations should be 
understood as proving that the exclusive DAR from Karmen [143, 144] and from LSND [145, 146, 147] 
are consistent with other measurements probing the nuclear matrix element. 

In fact the description of the GT transition between the ^^G and ground states (or to the 
other members of the T = 1 multiplett in the A = 12 nuclei including the ^^B ground state and the 
J = state in ^^G at 15.11 MeV) is a very demanding nuclear structure problem involving the subtle 
breaking of the SU(4) symmetry (in the limit of good SU(4) symmetry the GT transition vanishes). 
Glearly RPA correlations are insufficient and overestimate the transition strength by about a factor 
of 4 [150, 151, 126]. Hayes et al. performed shell model calculations in successively larger model 
spaces up to 6hw and gained two important insights into the problem [152]: i) enlarging the model 
space increasingly breaks the SU(4) symmetry. Relatedly the GT matrix elements increases moving 
towards the experimental value. However, this convergence process is quite slow, ii) The inclusion of 
a three-body interaction noticeable improves the agreement with data. The hndings of Ref. [152] have 
recently been conhrmed in no-core shell model (NGSM) calculations using chiral NN-I-3N interactions 
[153, 154]. In fact, the NSGM calculation of Maris et al. [154], performed in a model space with 8hw 
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Figure 5: Comparison of the exclusive ^‘^C{ue,e cross section data of the LSND collaboration 

[146] with shell model calculations performed with two different residual interactions (from [128]). 


now slightly overestimates the transition by about 15%. A good reproduction of the GT transition 
strength, and relatedly of the exclusive e“)^^Ngs has also been obtained in a diagonalization 

shell model calculation performed in a 3hw model space [128] using a novel residual interaction [155]. 
In fact, using a slightly quenched value for the axial-vector coupling constant Otsuka et al. [155] hnd 
very good agreement with the LSND exclusive cross sections (see Fig. 5). 

The inclusive ^^C(z/e, e“)^^N cross section induced by DAR neutrinos leads to many excited states in 
The typical momentum transfer in this process, g ~ 50 MeV, is noticeably larger than the average 
excitation energy, cj ~ 20 MeV. Under these conditions the cross section is dominated by forbidden 
transitions, but it is not too sensitive to details of the strength distribution. As discussed above, RPA 
studies should be adequate to reproduce the cross sections which is indeed the case as it is in details 
outlined in Ref. [18]. Further proof that cross sections, which are dominated by contributions from 
collective excitations, are well described by RPA calculations are given by studies of the muon capture 
on [150] and by inclusive electron scattering. The successful reproduction of muon capture rates by 
RPA calculations is not restricted to but has been extended to the data set across the full nuclear 
chart [156, 157]. 

The KARMEN collaboration has reported a measurement of the e~)^^Co cross section for 

DAR neutrinos, a = (2.56 ± 1.08 ± 0.43 ■ 10“'^*’) cm^ [148]. Table 1 compares this cross section with 
results from a hybrid model study and from different variants of the RPA. All calculations are consistent 
with the data within their relatively large uncertainty. We note, however, that the RPA cross sections 
show a spread of about a factor of 2.5. Looking at the results of the hybrid model in more details one 
hnds the largest contributions from the allowed transitions. The A = !■*■ part (1.12 x 10“^° cm^) has 
been calculated from the shell model GT strength distribution with the standard quenching factor and 
corrected for effects of hnite momentum transfer [149] . The A = 0^ contribution to the Isobaric Analog 
State at excitation energy of 3.5 MeV in ®®Go is readily calculated exploiting the appropriate Fermi sum 
rule yielding a cross section of 0.53 x 10“'^° cm^. Higher multipoles, which have been treated within an 
RPA approach, contribute a cross section of 0.73 x 10“^° cm^. The shell model A = I"*" cross section 
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Table 1: Comparison of ^®Fe(z/e,e )^®Co cross sections for DAR neutrinos calculated with different 
variants of the Random Phase Approximation, with the hybrid model and with data (exp) [148]. 


model cross section (in 10 cm^) 

QRPA [158] 

2.646 

PQRPA [158] 

1.973 

RPA [159] 

2.77 

QRPA [160] 

3.52 

RQRPA [119] 

1.4 

Hybrid [149] 

2.38 

exp [148] 

(2.56 ±1.08 ±0.43) 


is about 32% smaller than the one obtained within the RPA approach showing some sensitivity to the 
detailed GT distribution. 

Such GT distributions - as other relevant strength distributions for dipole transitions or Ml tran¬ 
sitions - serve as indirect constraints to validate the nuclear input into neutrino-nucleus cross section 
calculations. For ®®Fe Ref. [161] have measured the forward-angle cross section in a charge-exchange 
(p, n) reaction. These data are proportional to the GT_ distribution, which is of importance for the 
KARMEN measurement of the {ve-, e“) cross section. In Fig. 6 the measured forward-angle cross section 
is compared to the results obtained from a shell model calculation [29], where the latter has been folded 
with the experimental resolution. Overall a reasonable agreement is observed. The total GT_ strength 
has been calculated as 9.3 units, while Ref. [161] quotes 9.9 ± 2.4 units. In agreement with experi¬ 
ment, the calculation places the main strength at energies between 6 and 15 MeV, however, it misses 
some of the detailed fragmentation in the data. For low excitation energies the forward-angle cross 
sections have been recently determined using the (^He, t) charge-exchange reaction with a noticeably 
improved resolution compared to the pioneering (p, n) experiments [162, 163], again showing a rather 
good agreement with the shell model calculation. 

The left panel of Fig. 7 compares the GT_ strength calculated within a relativistic QRPA approach 
with a shell model calculation [164]. The shell model calculation uses a different residual interaction than 
adopted in [29], but yields very similar results. The total GT_ strength with the GXPFIJ interaction 
is 9.5 units, but it is also well reproduced by the QRPA study 9.7 units [164]. (Note that also the 
Ikeda sum rule is quenched.) The QRPA strength distribution has some strength at excitation energies 
below 5 MeV, similar to the shell model, and it places the largest portion of the strength in a resonance 
centred around 10 MeV. As expected, the shell model and experimental distributions show a signihcantly 
stronger fragmentation of this resonant strength which is not recovered in the QRPA study due to lack 
of correlations. 

As we have stressed above the detailed distribution of the GT_ strength matters in the description of 
neutrino-nucleus reactions if the incident neutrino energy is of the order of the nuclear excitation energy 
(considering the Q value of the reaction). This is nicely demonstrated in Fig. 7. For the distributions 
dominated by low energy neutrinos (at T = 2 and 3 MeV) the QRPA predicts slightly smaller cross 
sections than the hybrid model. But for the neutrino distributions with temperatures T > 4 MeV the 
cross sections are very similar indicating that the (average) neutrino energies are large enough that the 
cross section is only sensitive to the total strength and to the energy centroid of the distribution (which 
are quite similar in the QRPA and shell model), but not to the details of the distribution. Also the 
hybrid model calculation of Ref. [165], using a different shell model residual interaction and another 
RPA variant, yields cross sections which are very close to the one obtained in Ref. [164]. We note 
that the relative importance of correlations for the description of the total GT_ strength decreases with 
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Figure 6: Comparison of the GT strength distributions for ®®Fe obtained by {p,n) charge-exchange 
experiment [161] and a shell model calculation (folded with the experimental resolution). The data also 
include the Fermi contribution from the Isobaric Analog State (IAS) (from [29]). 
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Figure 7: (left panel) GT_ strength distribution for ^®Fe, calculated within a relativistic QRPA approach 
(RQRPA with interaction DD-ME2) and a shell model study with the GXPFIJ residual interaction. 
Both calculations have been folded with a Lorentzian of width 0.5 MeV. (right panel) Gomparison of 
(i^e; e“) cross sections on ^®Fe calculated within relativistic QRPA approaches with different quenchings 
of the axial-vector coupling constant (RHB-I-RQRPA) and within hybrid models where the GT_ strength 
has been calculated within the shell model and the other multipole contributions within the RPA. The 
cross sections have been folded with a Fermi-Dirac distribution with temperature T (corrected version 
from [164]). 


growing neutron excess as for nuclei with modest and large neutron excess this strength is dominated 
by the (quenched) 3(A^ — Z) part stemming from the Ikeda sum rule. 

We expect that the results found for ®®Fe applies quite generally: i.e. for charged-current {Uf,, e~) re¬ 
actions, RPA and hybrid model calculations will give reliable estimates for supernova neutrinos provided 
they give a reasonable account of the total GT_ strength and of the energy centroid. 

As the GT_|_ strength for nuclei with neutron excess is unconstrained by the Ikeda sum rule, its 
reproduction is a larger challenge to nuclear models than the GT_ strength. Importantly the shell 
model could be established as a reliable tool for such calculations in recent years [29, 31]. Using 
again ^®Fe as an example. Fig. 8 compares the experimental GT+ strength with various calculations 
performed within the shell model and the QRPA. At hrst we notice the signihcantly improved resolution 
obtained in the {t,^ He) charge-exchange experiment [168] (a similar resolution is achieved in the (d,^ He) 
experiments) compared to the pioneering work using the {n,p) reaction with an energy resolution of 
about 1 MeV [166, 167]. Secondly, the shell model gives a fair account of the total GT+ strength as 
well as of the distribution, although both residual interactions place slightly too much strength into 
the energy interval between 2 and 4 MeV. The shell model calculations, which agree quite well with 
each other, however, differ from the QRPA predictions. The latter yields noticeably more total strength 
than the data (about 5.5 units and 2.6 ± 0.8 units, respectively). Furthermore the lack of correlations 
in the QRPA has the consequence that not sufficient strength is pulled to excitation energies below 2 
MeV, in contrast to data and the shell model. Such shortcomings can lead to signihcant differences in 
the {i?e, e^) cross sections for low anti-neutrino energies. 

Satisfyingly the shell model has been proven to describe all measured GT_|_ distributions for pf shell 
nuclei quite well [169] . This ability makes it the method of choice to describe astrophysically important 
weak-interaction processes like electron capture or (z/g, e"*") reactions if the astrophysical conditions make 
the respective cross sections sensitive to a detailed description of the GT+ distribution. As an example. 
Fig. 9 compares electron captures rates calculated from the measured GT_|_ distributions for all pf 
shell nuclei, for which these data exist, and from QRPA and shell model strength distributions (for two 
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Figure 8: Comparison of GT+ strength for ®®Fe deduced from {n,p) [166, 167] and (t,^He) [168] 
charge-exchange experiments with distributions calculated within the QRPA approach (c) and by the 
shell model (d and e) with different residual interactions (KB3G and GXPFla). Panel f) shows the 
running sums of the strengths. The QRPA results have been divided by a factor of 3 (from [168]). 
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Figure 9: Comparison of electron-capture rates for 13 p/-shell nuclei calculated on the basis of theoretical 
(connected by lines) and experimental (indicated by markers) GT strength distributions. All rates are 
plotted relative to those calculated by the shell-model with the KB3G interaction. The rates are shown 
for two combinations of stellar density and temperature: (a) pY^. = 10^ g cm“^ and T = 3 x 10® K and 
(b) pYf. = 10® g cm“^ and T = 10 x 10® K (from [169]). 


different residual interactions). The astrophysical conditions in the right panel correspond to the early 
onset of collapse where the electron Fermi energy is of the same order of magnitude as the reaction 
Q value making the capture rate sensitive to details of the GT_|_ distribution. Fig. 9 shows a very 
satisfying agreement between the rates obtained from data and shell model distributions. The capture 
rates obtained from the QRPA calculations agrees usually also within a factor of 2. 

The sensitivity of the rates to the details of GT+ distributions decreases with increasing Fermi 
energy (or density). Fig. 9 conhrms this fact by showing a noticeably larger variation between the 
different calculations and the data for supernova conditions at smaller density. In particular two points 
are noteworthy: i) the better energy resolution achieved in the (d,^ He) and (f,^ He) reactions usually 
leads to smaller rates than those obtained from the (u,p) data, in good agreement with the shell model; 
ii) the missing fragmentation of the QRPA strengths can lead to large differences in either directions. 
We stress, however, that the pf shell nuclei shown in the figure have small abundances at the relatively 
low densities and hence do not contribute much to the overall electron capture rate. 

A similar study as performed for electron capture in Ref. [169] would be also desirable for (z/g, e"*") 
cross sections for supernova neutrinos. As energies of supernova Ve neutrinos are larger than the electron 
Fermi energies for both collapse situations depicted in Fig. 9. we expect from the general considerations 
above that the respective (z/g, e"*") cross sections are less sensitive to the details of the GT_|_ distribution 
than shown at the conditions of Fig. 9 and that the shell model, and also the QRPA, should give a fair 
account of the GT contribution to the cross sections, supplemented by higher multipole contributions 
obtained from RPA calculations. This statement is likely correct for the chosen pf shell nuclei, but 
it can unfortunately not be generalized. For these pf shell nuclei the neutron number is smaller than 
40 and hence GT+ transitions, in which a proton is changed into a neutron, is even possible within 
the IPM. However, for nuclei with proton numbers below a major oscillator shell closure (say Z < 20 
or Z < 40) and neutron numbers above {N > 2Q or N > 40, respectively) GT_|_ transitions are Pauli 
blocked in the IPM. However, they can open up by correlations which move nucleons across the shell 
closure. Thus such correlations can enable GT+ transitions for protons in the upper oscillator shell 
or into neutron holes in the lower oscillator shell. Detailed calculations have shown that describing 
cross-shell correlations is a slowly converging process requiring many-particle many-hole configurations 
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Figure 10: (left panel) Running GT_|_ sum for different values of the particle-hole excitations across 
the = 40 shell gap. (right panel) Electron capture rates calculated with normalized GT strength 
distributions at different truncation levels (from [40]). 


usually beyond those considered in QRPA studies [170, 171]. To illustrate this point we use the electron 
capture rate on ^®Se with = 42 and Z = 34 as an example (see Section 2). the left panel of Fig. 10 
shows the shell model results for the ^®Se GT+ strength. Note that the shell model space is adequate 
to describe the correlations of valence nucleons across the shell gap, it lacks, however, some spin-orbit 
partners making the use of the universal quenching factor inappropriate. In Ref. [40] the quenching 
has been simply chosen to match the total experimental GT_|_ strength [38]. This choice does not affect 
our argument here which is the evolution of the strength with truncation level; i.e. with the number of 
nucleons allowed to be promoted from the pf shell across the shell gap by correlations. In particular, 
truncation level t = 2 includes up to two-particle two-hole correlations like in QRPA calculations. We 
observe from Fig. 10 that, with increasing correlations, i) the total GT+ strength grows, mainly due to 
an increase of neutron holes in the pf shell, and ii) more strength is shifted to lower excitation energies. 
This has an important effect on the capture rates. This is illustrated in the right panel of Fig. 10 at 
conditions with electron Fermi energies of nearly 10 MeV, where the rates differ by up to an order 
of magnitude between the calculation with truncation level t = 2 and the converged calculation. We 
mention that in the IPM the GT_|_ strength and the rates vanish. As the average energies of supernova 
z/g neutrinos are slightly larger than the electron Fermi energy in our example, we expect that the 
differences in [ug, e"*") cross sections for supernova neutrinos will be still appreciable, but somewhat 
smaller. 

The only neutral-current cross section measured for nuclei beyond the deuteron is the transition from 
the ground state to the T = 1 state at 15.11 MeV in ^^G [144, 172]. This is the Isobaric Analog State of 
the and ground states. Thus from a nuclear structure point of view this measurement yields the 
same informations as the exclusive charged-current reactions on ^^G discussed above. Hence we have to 
use indirect experimental information to validate cross section calculationsi for inelastic neutrino-nucleus 
scattering. Due to the relatively low energies of supernova neutrinos allowed transitions dominate 
inelastic (z/, z/') cross sections and the relevant nuclear quantity is the GTq strength distribution. This 
strength can be experimentally studied for example in {p,p') experiments. However it is also intimitely 
related to Ml excitations of nuclei which have been studied with great precision by inelastic electron 
scattering [173]. We will now use such electron scattering data to show that the shell model is capable 
to describe the allowed response needed to determine inelastic cross sections for supernova neutrinos. 
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Ml transitions are mediated by the operator 


0(M1) 



'Yy9i{k)l{k) + gs{k)s{k)]iiN 

k 


( 19 ) 


where the sum runs over all nucleons. The orbital and spin gyromagnetic factors are given by 5 ^/ = 1, 
gs = 5.585694675(57) for protons and gi = 0, gs = —3.82608545(90) for neutrons [174]; /iat is the nuclear 
magneton. Using isospin quantum numbers = T 3/2 = ±1/2 for protons and neutrons, respectively, 
Eq. (19) can be rewritten in isovector and isoscalar parts. Due to a strong cancellation of the gf-factors 
in the isoscalar part, the Ml transitions are of dominant isovector nature. The respective isovector Ml 
operator is given by 


0(Afl)„ = J^ J2l^{k)Uk) + (9? - (20) 

k 

We note that the spin part of the isovector Ml operator is the zero component of the GT operator, 

0{GTo) = Y, Mk)Uk) = Y (21) 

k k 

however, enhanced by the factor ^y3J^{g^ — 5'/)/2 = 2.2993. The allowed contribution to the in¬ 
elastic neutrino-nucleus scattering at low energies, where hnite momentum transfer corrections can be 
neglected, from an initial nuclear state {i) to a hnal state (/) is then given by 

* ^ /) = {/II <T(A.)t.(m|i> 1^ (22) 

7r(2 ±1) ^ 

where is the energy of the scattered neutrino and u is the difference between final and initial nuclear 
energies. The nuclear dependence is contained in the B{GTq) = (y'^|(/|| '^k^(k)tz{k)\\i)\'^/{2Ji + 1) 
reduced transition probability between the initial and final nuclear states. In shell model calculations of 
the Ml strength the spin part is quenched by the same factor used for studies of the GT strengths [175]. 

Thus, experimental Ml data yield the desired GTq information, required to determine inelastic 
neutrino scattering on nuclei at supernova energies, to the extent that the isoscalar and orbital pieces 
present in the M1 operator can be neglected. On general grounds one expects, as stated above, that the 
isovector component dominates over the isoscalar piece. Furthermore, it is wellknown that the major 
strength of the orbital and spin Ml responses are energetically well separated in the nucleus. In pf shell 
nuclei, which are of interest for supernova neutrino-nucleus scattering, the orbital strength is located at 
excitation energies E* ~ 2 — 4 MeV [176], while the spin Ml strength is concentrated between 7 and 
11 MeV. A separation of spin and orbital pieces is further facilitated by the fact that the orbital part 
is strongly related to nuclear deformation [177]. For example, the famous scissors mode [178], which 
is the collective orbital Ml excitation, has been detected in well-deformed nuclei like ®®Fe [179]. Thus 
one can expect that in spherical nuclei the orbital Ml response is not only energetically well separated 
from the spin part, but also strongly suppressed. 

Examples of spherical pf shell nuclei are ^°Ti, ®^Cr and ®^Fe. In Fig. 11 we show the total Ml 
strength and its spin and orbital contributions from a shell model calculation for ^^Cr and compare it 
to the experimental data [180] . We observe a very good agreement between data and calculation for the 
Ml strength distribution, with the strength mainly located at excitation energies between 7-10 MeV. We 
also find that the calculated Ml strength is dominated by its spin contribution which exceeds the orbital 
part by more than an order of magnitude. Finally, the calculation yields no indication for a scissors mode 
or any orbital strength at low excitation energies, as expected for spherical nuclei. Thus, the calculation 
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Figure 11: Comparison of experimental ®^Cr Ml strength distribution (bottom panel) with shell-model 
result (upper panels) showing the orbital, spin components and the total. Note the different scales of 
the ordinate for the spin and orbital pieces, respectively. The shell-model calculations show isospin 
T = 2 states as blue spikes and T = 3 as red spikes (adapted from [180]). 
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Figure 12: Neutrino-nucleus cross sections on ®^Cr calculated from the Ml data (black line) and the 
shell-model GTq distribution (red line). The blue line shows the cross section for a temperature of 
0.8 MeV (adapted from [180]). 


strongly supports the argument of Ref. [180] that, except for an overall normalization, the Ml strength 
for spherical nuclei reflects in a good approximation the GTq distribution. Based on this assumption. 
Ref. [180] has calculated the allowed contributions to the inelastic neutrino scattering on ®^Gr from both 
the Ml data (properly normalized to represent the GTq strength) and the GTq strength calculated in 
the shell model. Both cross sections are compared in Fig. 12 and show excellent agreement for scattering 
from the ground state. For the cross section at finite temperature contributions of thermally populated 
states have to be included. These have been calculated from appropriate shell model GTq strength 
distributions and added to the ground state cross section as described above, eq. (15). Fig. 12 indicates 
that the contributions from excited states have noticeable effects at low neutrino energies. This is 
particularly enhanced for spherical nuclei like ®^Gr where there is no significant GTq strength at low 
nuclear excitation energies. For neutrino energies which are sufficient to excite the GTq resonance, 
related to the spin part in the Ml strength, from the ground state the finite temperature corrections 
become negligible. For applications to supernova neutrinos contributions from forbidden transitions 
have to be added to the cross section. These can be calculated within RPA approaches and become 
increasingly relevant with growing neutrino energies. 

The importance of a detailed reproduction can be deduced from a study presented in Ref. [134]. 
These authors performed a comparison of inelastic neutrino-nucleus cross sections for ^®Fe calculated 
within the hybrid model (shell model for allowed contributions and RPA for other multipoles, [165] and 
with studies in which all multipole contributions are determined within different variants of the QRPA 
[181, 182, 134]. The cross sections are summarized in Table 2 for different incident neutrino energies. 
We note some differences in the cross section predictions between the QRPA calculations. As mentioned 
in Ref. [182] these are related to the choices of single particle energies and the residual interaction, and 
should become less important with increasing neutrino energy. However, the calculations also treated 
the quenching of the various multipole contributions differently: only quenching the GTq contribution 
(like in the hybrid model) [134] or quenching all multipole contributions [182]. As Ref. [134] took 
special care in reproducing the centroid and strength of the GTq distribution and of the dipole giant 
resonance by slight adjustment of the residual interaction, this study allows a detailed comparison to 
the hybrid model calculation [165]. We find indeed a very good agreement between the cross sections 
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except at the lowest neutrino energy = 10 MeV. Here the cross section is sensitive to a detailed 
description of the GTq strength; the shell model describes the fragmentation of the strength while the 
QRPA calculations predict it to be quite strongly concentrated around excitation energies of = 10 
MeV (see Figs. 15 and 13 below). 


Table 2: Comparison of cross sections (in units of 10“^^ cm^) for inelastic neutrino scattering on the 
ground state of ^®Fe calculated within the hybrid model [165] and different variants of the QRPA. The 
incident neutrino energy is Ey. Exponents are given in parenthesis (from [134]). 


Ey [MeV] 

hybrid model [165] 

QRPA [134] 

QRPA [182] 

QRPA [181] 

10 

1.91(-1) 

1.69(-2) 

1.87(-1) 

1 .01(+0) 

20 

6.90(+0) 

5.64(+0) 

9.78(+0) 

5.79(+0) 

30 

2.85(+l) 

2.41(+1) 

4.08(+l) 

1.87(+1) 

40 

7.86(+l) 

6.65(+l) 

1.05(+2) 

5.51(+1) 

50 

1.72(+2) 

1.49(+2) 

2.16(+2) 

1.43(+2) 

60 

3.20(+2) 

2.87(+2) 

3.89(+2) 

3.09(+2) 

70 

5.25(+2) 

4.83(+2) 

6.33(+2) 

5.63(+2) 

80 

7.89(+2) 

7.36(+2) 

9.59(+2) 

8.82(+2) 

90 

l.ll(+3) 

1.03(+3) 

1.38(+3) 

1.22(+3) 

100 

1.49(+3) 

1.36(+3) 

1.92(+3) 

1.52(+3) 


4 Applications to supernova dynamics 

As we have stressed before neutrino-induced processes play a crucial role in the dynamics of core¬ 
collapse supernovae. In accordance with the main theme of this review we will focus here on two 
types of neutrino-nucleus reactions which both have recently for the hrst time been incorporated in 
supernova simulations and their impact on the dynamics has been explored. The hrst subsection deals 
with inelastic neutrino-nucleus scattering. We have outlined in the last section how progress in nuclear 
modelling, constrained by experimental data from inelastic electron scattering and charge-exchange 
experiments, has made possible to reliably estimate respective cross sections for scattering on ground 
states, despite the absence of direct data, and has shown ways to extend the calculations to hnite 
temperature as is appropriate for supernova simulations. In the hrst subsection we will discuss the 
results of the nuclear model calculations, focussing on the ehects introduced by hnite temperature, and 
then turn to the impact elastic neutrino-nucleus scattering has in supernova simulations. In the second 
subsection we deal with nuclear deexcitation by neutrino pairs; i.e. the transition of a nucleus from a 
state at higher energy to a state at lower energy by neutrino pair production. Obviously this is only 
possible at hnite temperature and can therefore not be directly studied in the laboratory. However, 
from a nuclear structure point of view this neutral current process is dominated by the same allowed 
Gamow-Teller and forbidden response as inelastic neutrino scattering, but in diherent kinematics. We 
will discuss two quite distinct approaches to derive the respective cross sections. Both yield, however, 
similar results. Supernova simulations, which include nuclear deexcitation by neutrino pair emission, 
indicate that this process does not impact the dynamics of the collapse. The process is, however, a 
source of and z/g neutrinos during the collapse phase. 
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4.1 Inelastic neutrino-nucleus scattering and its impact in supernova sim¬ 
ulations 

Inelastic neutrino-nucleus scattering, Ai -\- v ^ Af -\- , changes the nucleus A from an initial state Ai 

with energy Ei to a hnal state Af with energy Ef. In this process the neutrino either gains energy if 
Ef < Ei OT it looses energy if Ef > E^. The energy exchange between the nucleus and the neutrino 
can have potentially interesting consequences for the supernova dynamics. Inelastic neutrino-nucleus 
scattering might i) transfer energy to nuclei outside the shock front after bounce before arrival of the 
shock [183] enabling an easier revival of the shock, ii) add to inelastic neutrino-electron scattering 
speeding up the thermalization of neutrinos with matter after neutrino trapping during the collapse 
phase, and hi) change the neutrino opacity which will in turn modify the spectra of neutrinos released 
in the supernova and subsequently observed by earhbound detectors. 

Investigations of inelastic neutrino-nucleus scattering in SN simulations has been pioneered in an 
exploratory study [184], approximating the matter composition by a representative nucleus, ®®Fe. The 
reaction cross sections were based on a nuclear model appropriate for temperature T = 0, combining 
a truncated shell model evaluation of the allowed Gamow-Teller response to the cross section with 
estimates of forbidden components derived from the Goldhaber-Teller model. The study concluded 
that inelastic neutrino-nucleus scattering rates can compete with those of neutrino-electron scattering 
at moderate and high neutrino energies {E^ > 25MeV), while they are significantly smaller for low E^. 
No significant effects of neutrino-nucleus interactions on the stalled shock by preheating the accreted 
matter were found in [184]. 

Approximating the composition of supernova matter by the ground state of the even-even nucleus 
®®Fe is too simple an assumption for the calculation of inelastic neutrino-nucleus interaction rates for 
two reasons, i) Gamow-Teller transitions, which as we recall dominate the cross sections at the neutrino 
energies of interest, can only connect the O’*" ground state of ^®Fe to hnal states with spin and parity 
J = I"*". As the lowest I"*" state in ^®Fe is at an excitation energy oi E^ = ?> MeV, there exists a threshold 
for inelastic neutrino scattering on ®®Fe (similarly on other even-even nuclei) and the cross sections are 
rather small for low neutrino energies. Supernova matter consists of a mixture of many nuclei with even 
and odd proton and neutron numbers. Since odd-A and odd-odd (i.e., with odd proton and neutron 
numbers) nuclei miss the strong pairing gap that lowers the ground state in even-even nuclei relative 
to excited states, and have usually J* 7 ^ 0, Gamow-Teller transitions from the ground state to levels at 
rather low excitation energies are possible, reducing the threshold for inelastic neutrino scattering on 
the ground state and generally increasing the cross sections at low E^. ii) More importantly, however, 
supernova matter has a non-zero temperature of order 1 MeV or higher, requiring the description of 
nuclei as a thermal ensemble. This completely removes the energy threshold for inelastic neutrino 
scattering, as we have shown in Section 3. 

The hnite temperature effect on the inelastic neutrino cross section off the two even-even nuclei 
®®Fe and ®^Ge has been in detailed explored recently by Dzhioev and collaborators [134] using the 
TQRPA approach. As their general hndings for the two nuclei are similar, we concentrate on ^®Fe 
here. The temperature evolution of the GTq strength is exhibited in Fig. 13. The ground state 
(T = 0) distribution shows that the strength is overwhelmingly concentrated at excitation energies 
around 10 MeV, corresponding to a spin excitation. The sttrength at low energies around 4 MeV 
is signihcantly lower. We note that these states are more pronouncedly seen in Ml data as they 
predominantly reflect orbital excitations, including the scissors mode at 3.45 MeV. At hnite temperature 
the GTq distributions show several effects [134] which are relevant for the subsequent discussion of 
neutrino-nucleus scattering. For positive neutrino energy transfer (i? > 0 in Fig. 13) the main effect 
is a shift of the strength towards smaller energies, by about 1.5 MeV. This shift is attributed to the 
pairing phase transition [134] : at temperatures above the critical one no extra energy has to be invested 
into breaking Gooper pairs and as a consequence the GTq strength moves to lower energies. This 
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Figure 13: Temperature evolution of the GTq strength distributions for ®®Fe (left panels) and ®^Ge (right 
panels) calculated in the TQRPA approach. The x axis denotes the neutrino energy transfer to the 
nucleus. The red arrows indicate the threshold for inelastic neutrino scattering on ®®Fe at temperature 
T = 0 which is given by the excitation energy of the lowest J = state. Here the experimental value 
3.12 MeV is chosen, which is slightly lower than the calculated value (from [134]). 
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Figure 14: (upper panels): Cross sections for inelastic neutrino-nucleus scattering on ®®Fe (left) and ®^Ge 
(right) at different temperatures: T = 0 (full line), T = 0.86 MeV (blue dashed line), T = 1.29 MeV 
(brown dot-dashed line), T = 1.72 MeV (red double-dot-dashed line), (lower panels): Contribution of 
the forbidden multipoles to the cross sections at selected temperatures (from [134]). 


behavior will not be present in calculations which are based on Brink’s hypothesis. Furthermore, at 
hnite temperature strength appears now also at negative neutrino energy transfers {E < 0). This 
relates to thermally populated nuclear states which decay to states at lower energies. Obviously the 
downward transitions increase noticeably with temperature. The strength observed a.t E = —9 MeV can 
be attributed to the deexcitation of the GTq resonance. Dzhioev et al. point out that the assumption 
of Brink’s hypothesis also affects the downward strength distribution. In their approach the shift of the 
upward strengths to lower energies due to the melting of pairing makes it easier to thermally populate 
these states which, in turn, leads to an enhanced downward strength as well. Furthermore we remark 
that the smearing of the Fermi surface at hnite temperature allows for a larger fragmentation of the 
strength, which affects more the low lying transitions (and ®^Ge). 

The Gamow-Teller transitions and contributions from forbidden multipoles up to A = 3 have been 
used to calculate inelastic neutrino-nucleus cross sections for ^®Fe and ®^Ge at different temperatures 
[134]. The results are shown in Fig. 14. There are several interesting observations, i) At T = 0, i.e. 
using the ground state transitions, the cross section drops at low neutrino energies. This is related to 
the gap in the GT distribution requiring neutrinos to have a certain minimum energy (about 4 MeV 
in the TQRPA ®®Fe calculation) to excite the lowest J = state, ii) There is no gap in the cross 
sections at the hnite temperatures mainly caused by the fact that deexciation of thermally populated 
nuclear states contributes to the cross section at all neutrino energies. Obviously the contributions of 
the thermally populated states increase with temperature, so does the scattering cross section at low 
energies, iii) With increasing neutrino energies the cross sections at diherent temperatures converge. 
At neutrino energies E < 15 MeV, excitation of the GTq resonance dominates the cross section. At 
the energies shown in Fig. 14 the shift of the centroid energy due to the melting of pairing leads to 
a small increase in cross section at the hnite temperatures compared to the T = 0 case. Due to the 
energy dependence of the phase space factor the relative enhancement of the cross section caused by 
this energy shift of the centroid decreases with increasing neutrino energy, iv) At the neutrino energies 
shown, which are typical for supernova neutrinos, the cross sections is dominated by the GT contribution 
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Figure 15: Gaussian-smoothed neutral Gamow-Teller (GTq) distributions from the ground state in Fe 
isotopes. Solid lines indicate AT = 0 strength, and dashed lines AT = 1. is the nuclear excitation 
energy (from [130]). 


at all temperatures. Forbidden transitions contribute less than 10%. 

We mentioned above briefly that the GT distributions and the inelastic neutrino-nucleus cross sec¬ 
tions at supernova neutrino energies are quite similar for ®®Fe and ®^Ge. This is actually more globally 
true. Neutral-current (z/, u') cross sections, folded over appropriate supernova neutrino energy distribu¬ 
tions and calculated for T = 0, show relatively small variations from nucleus to nucleus [9, 104, 105]. 
The reason is related to the fact that the cross section at the supernova energies is dominated by GT 
transitions and that the GT strength is strongly concentrated in a resonance at excitation energies 
around 9-10 MeV for medium-mass and heavy nuclei. To support this point we show the GTq distri¬ 
butions from shell model calculations [130] for the iron isotopes 52-6ipg pjg_ [185]). In 

fact for all isotopes the dominant part of the strength resides in a resonance around 10 MeV. Also the 
total resonant strength does not vary much. Very similar strength functions are found for the mangan, 
cobalt and nickel isotopic chains [130]. It is worth mentioning that the GTq distributions for this set of 
Mn, Fe, Go and Ni isotopes does not depend on the pairing structure of the ground state (even-even, 
odd-A, odd-odd). 
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We recall that the hybrid model (introduced in Section 3) explicitly incorporates Brink’s hypothesis 
for the downscattering part of the cross section. Hence it does not account for the energy shift in the 
GTo strength due to the melting of pairing, in contrast to the TQRPA approach. As a consequence 
the TQRPA cross sections for ^®Fe (the only nucleus which has been studied in both approaches) are 
somewhat larger (by about a factor of 2-3) at < 10 MeV than those calculated in the hybrid model 
[134]. For higher energies the cross sections of the two models become nearly identical. 

Juoadagalvis et al. have calculated inelastic neutrino scattering cross sections on Mn, Fe, Co and Ni 
isotopes at hnite temperatures [130]. The study has been performed within the hybrid model (taking 
the GTq distributions from the shell model and other multipole contributions from RPA calculations). 
The hnite temperature effects have been derived by inverting the shell model GTq distributions for the 
lowest excited states. The cross sections for the iron isotopes are shown in Fig. 16. We again observe the 
same trend in the cross sections as discussed above for TQRPA calculations. At low neutrino energies 
the cross sections is given by the contributions arizing from deexcitation of thermally populated states 
and there are small variations in the cross section between the different isotopes. Once the neutrino 
energy is large enough to excite the GTq resonance the cross sections for all isotopes become very similar 
(which is due to the similarity of the GTq distributions) and hnite-temperature effects become rather 
small. At energies < 20 MeV forbidden contributions are at least an order of magnitude smaller 
than the GT contribution, become, however, increasingly more important with energy. 

The effect of inelastic neutrino-nucleus scattering on the supernova dynamics has been in details 
explored in Ref. [5] by simulating the collapse of a 15 Mq progenitor star (model sl5a28 of [186]). The 
simulation has been performed in spherical symmetry with the neutrino-hydrodynamics code VERTEX 
of the Garching group [89, 90] using a state-of-the-art description of the interactions of neutrinos and 
antineutrinos of all flavors as dehned in [90, 35, 84]. 

The inelastic neutrino-nucleus scattering cross sections were derived on the basis of the hybrid model 
for a pool of about 50 nuclei covering the Z = 24 — 28 chains [130, 5]. Importantly, for temperatures 
T > 1 MeV, these calculations do not show large variations between the individual cross sections, for 
the reasons explained above. This justihes to approximate the relevant neutrino-nucleus cross sections 
by replacing the supernova matter composition by the results drawn from the restricted pool of nuclei. 
Hence, Ref. [5] dehned an average neutrino-nucleus cross section as V, where the 

sum is restricted to the pool of nuclei for which individual cross sections exist. The abundances of 
the individual nuclei V were determined from an appropriate NSE matter distribution [187]. A rate 
table of the cross section has been derived as a function of temperature, density and electron-to-baryon 
ratio Ye allowing for cross section determination at all encountered astrophysical conditions during the 
simulation. As an example the pool-averaged total cross section is shown in Fig. 17 at 3 different 
temperatures. (The densities and Ye values are different at the different temperatures. They inhuence 
the average cross section by change of the NSE distribution with very minor impact.) The average cross 
section rehects the general trends discussed at length above. 

Keeping track of neutrino energies is a crucial ingredient in supernova simulations. Inelastic neutrino 
scattering on nuclei changes the energy of the incident neutrino. This change depends on the incident 
energy and on temperature. The right panel in Fig. 17 shows the normalized outgoing neutrino spectra 
at different initial neutrino energies and for temperature T = 1 MeV. For E^ = Q and 10 MeV down¬ 
scattering contributes signihcantly to the cross section; i.e., the de-excitation of thermally populated 
nuclear levels produces neutrinos with E'^ > 10 MeV. Down-scattering becomes essentially irrelevant at 
higher neutrino energies. For Ey = 30 MeV the cross section is dominated by the excitation of the GTq 
centroid, giving rise to a cross section peak around = 20 MeV. For even higher neutrino energies 
forbidden transitions contribute noticeably to the cross section. The peaks in the differential cross 
sections for E^ = 50 MeV correspond to excitations of the centroids of the GTq and dipole transition 
strengths. With increasing population of excited nuclear states the down-scattering contributions get 
relatively enhanced with increasing temperature [130, 134]. 
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Figure 16: Neutral current neutrino-nucleus inelastic cross-sections for Fe isotopes. The dashed lines 
show shell model contributions at three different temperatures (T = 0.86 MeV dashed line, 1.29 MeV 
dash-dotted line, and 1.72 MeV dash double-dotted line); the dotted line is the RPA contribution; 
the solid line corresponds to the total cross sections. At low neutrino energies the total cross section 
coincides with the SM contributions (from [130]). 
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Figure 17: (left panel) Inelastic neutrino scattering cross sections for supernova conditions (see text) 
at 3 different temperatures, (right panels) Normalized neutrino spectra derived from inelastic neutrino 
scattering on a supernova matter composition for 4 different initial neutrino energies and at temperature 
T = 1 MeV. For a comparison the dotted lines in the various panels give the respective quantities for 
the individual nucleus ^®Fe at T = 0 (based on the simulations presented in [5]). 



Figure 18: Normalized Ug number spectra radiated during the shock breakout burst as seen by a distant 
observer at rest. Results are shown for simulations with three different nuclear EoSs. For better 
comparison of the strongly time-dependent spectra during this evolution phase, integration in a window 
of 8 ms around the peak luminosity was performed. Inelastic neutrino scattering off nuclei (dashed 
lines) leads mostly to energy losses of high-energy neutrinos and thus reduces the high-energy tails of 
the spectra. The vertical line marks the mean spectral energy (from [5]). 


36 






















Table 3: Electron neutrino cross sections for scattering off electrons and for charged-current interactions 
with nuclei of different detector materials, averaged over the supernova neutrino spectrum of Fig. 18, 
using the results obtained with the Lattimer-Swesty EoS. The different cross sections were taken from: 
e [188], d [189], [151], [190], ^^Ar [18], ^^Fe and ^ospb [104] (from [5]). 


Material 

(a) (10 
With IS 

cm^) 

Without IS 

Reduction 

e 

0.114 

0.123 

7% 

d 

5.92 

7.08 

16% 

12C 

0.098 

0.17 

43% 

(N,3) 

0.089 

0.15 

41% 

160 

0.013 

0.031 

58% 

40 Ar 

17.1 

21.5 

20 % 

56Fe 

8.81 

12.0 

27% 

208pb 

147.2 

201.2 

27% 


Despite causing a higher neutrino opacity by increasing the neutrino-matter interaction, inelastic 
neutrino-nucleus scattering has virtually no impact on the supernova evolution [5]. As suggested by 
Haxton [183] the scattering increases the neutrino-matter coupling and thus the total energy transfer 
rate mainly to matter (nuclei and electrons) ahead of the shock by a factor 2-3 during the hrst 50 
ms after shock bounce. However, the time period during which this additional heating acts on the 
supersonically infalling matter is too short to have consequences on the supernova simulation. The 
simulation of [5] thus conhrms the expectation of Bruenn and Haxton [184] that neutrino preheating of 
matter ahead of the shock does not add signihcantly to the shock revival. 

Inelastic neutrino-nucleus scattering has, however, a signihcant influence on the spectra of neutrinos 
emitted after bounce, in particularly during the period of the so-called z/g neutrinoburst (which, as 
we recall, is caused by fast electron captures on free protons after photodissociation of nuclei by the 
shock, see Section 2). During this period matter ahead of the shock still consists of heavy nuclei in 
NSE. Furthermore the density is still high with a relatively large optical depth for neutrinos. Under 
these conditions neutrinos interact frequently with nuclei by inelastic scattering. Due to the energy 
dependence of the cross section the scattering occurs more often for high-energy neutrinos, which will 
spend part of their energy by mainly exciting the GTq or giant dipole resonances (see Fig. 17). In 
fact, inelastic neutrino-nucleus scattering shifts the mean spectral energy of the v^. burst spectrum only 
slightly by about 0.5 MeV, but it strongly suppresses the high-energy tail of the spectrum. This is 
shown in Fig. 18 which compares the z/g burst spectrum in supernova simulations with and without 
consideration of inelastic neutrino-nucleus scattering. A similar suppression of the high-energy tail, 
caused by inelastic neutrino-nucleus scattering, is also observed for z7g and neutrinos during this 
early phase of the explosion. At later times inelastic neutrino scattering on nuclei becomes negligible 
for the neutrino emission spectra as the abundance of heavy nuclei has dropped drastically due to 
dissociation of nuclei 

The observation of the z/g burst neutrinos is one of the aims of supernova neutrino detectors [191, 
192, 193]. The change of the z/g burst spectrum due to inelastic neutrino-nucleus scattering will have 
consequences on the expected event rates for such detectors. We note again that, due to the energy 
dependence of neutrino-nucleus cross sections, the weight of high-energy neutrinos gets enhanced in 
the event rates. In Table 3 we compare the relevant detection cross sections for several target nuclei 
calculated for the v^. burst spectra obtained with and without consideration of inelastic neutrino-nucleus 
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scattering. The largest reduction occurs for and which are detector materials in Borexino, 
MiniBooNe, KamLAND, SNO+, and Super-Kamiokande. For these nuclei only neutrinos with relatively 
high energies {E > 17 MeV for and E > lb MeV for can trigger charged-current reactions. 
As the neutrino scattering cross section on electrons increases linearly with energy, the reduction of 
the event rate for electrons, is less. Ref [5] estimates the improved event rates (including the effect 
of inelastic neutrino-nucleus scattering) arizing from the z/g burst of a supernova going off at 10 kpc 
distance. 

In summary, supernova simulations imply that inelastic neutrino-nucleus scattering has no impact 
on the supernova evolution, but noticeably changes the Vf. burst spectrum. The latter effect is mainly 
caused by scattering on nuclei induced by high energy neutrinos, for which the respective cross sections 
are likely under control and can be calculated with sufficient accuracy, also at hnite temperature. 
Finite temperature corrections signihcantly affect the cross sections for low-energy neutrinos. Here 
recent developments achieved within the TQRPA approach [134] are very promising and show that 
the cross section estimates based on the hnite-temperature treatment in the hybrid model [130] are 
probably somewhat too small. This uncertainty in the low-energy neutrino-nucleus cross section has, 
however, probably no impact on the supernova evolution or the emitted neutrino spectra. 

4.2 Deexcitation of hot nuclei by neutrino-pair emission in supernova en¬ 
vironment 

During the hnal phase of the collapse the temperature in the inner core can reach values above 1.5 MeV. 
Under such conditions the average excitation energy of a nucleus with mass number A ~ 120, (a repre¬ 
sentative value for the matter composition of heavy nuclei) is about {E^) = 35 MeV. It has been argued 
that the deexcitation of such highly excited states can be a signihcant neutrino-pair producing process 
[194] which can compete with electron-positron annihilation or nucleon-nucleon bremsstrahlung. As the 
presence of electron neutrinos does not block the production of muon and tau neutrino-antineutrino 
pairs, the nuclear deexcitation process might further reduce the entropy of the collapsing core if the 
neutrinos produced by the process can leave the core during the dynamical timescale of the collapse 
[194] . Recently nuclear deexcitation by neutrino-pair emission has, for the hrst time, been incorporated 
into a supernova simulation [195] . We will here briefly summarize the results of this study. 

The formalism how to derive the rate for nuclear dexxecitation, Ai ^ Af + v + v, and its inverse 
reaction neutrino-pair absorption, Ai + v+V —)■ Af, and how to incorporate it into supernova simulations 
is presented in Ref. [195]. For our discussion here it is sufficient to recall that both processes are 
mediated by neutral current. In fact for neutrino-pair absorption the nuclear input is given by the 
same allowed and forbidden strength functions which determine inelastic neutrino-nucleus scattering. 
Our discussions above have shown that the hybrid model and the TQRPA calculations both indicate 
that these allowed and forbidden strengths are both concentrated around the GTq and giant dipole 
resonances. Furthermore the calculations (and the available Ml data) show that the position and 
total strengths do only mildly vary between nuclei. The authors of Ref. [195] used these hndings to 
approximate the strength for the absorption process by: 

^‘^"'^(A) = SAQAiA.) + SpgFiA.) (23) 

where A. = Ey + E^ is the energy sum of the neutrino-antineutrino pair. The normalized strength 
functions qa for allowed transitions and gp for forbidden transitions have been chosen as Gaussians. 
The centroids of these Gaussians ha, have been chosen guided by the shell model and TQRPA 
calculations for allowed transitions (/ia = 9 MeV, [134, 130, 185]) and by RPA calculations for forbidden 
transitions {gp = 22 MeV, [165, 134]). Assuming that the strengths for excited states is somewhat more 
fragmented and wider spread in energy the standard deviations ua, cp have been adopted slightly larger 
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Figure 19: Strength for absorption (top panels) and emission (bottom panels) of neutrino pairs based 
on the Fermi-gas model of ref. [194] (green lines) and the Gaussian ansatz, eq. (23), (red lines) at 
temperatures of T = 0.7 MeV (left panel), T = 1.5 MeV (middle panel), and T = 3.0 MeV (right 
panel). The average thermal energy at the different temperatures is given above the lower panels (from 

|195]). 


than for the ground state; i.e = 5 MeV, ap = 7 MeV. Finally the total allowed strength, S'^ = 5 and 
the total forbidden strength, Sp = 7, has been adopted from shell model [130, 185] and RPA calculations 
[165]. Finally we mention that the approximate ansatz for the absorption strength function, eq. (23), 
explicitly assumes the Brink hypothesis which makes the strength independent of the initial nuclear 
excitation energy and consequently the absorption strength becomes independent of temperature. 

The process of relevance for supernova simulations is nuclear deexcitation by neutrino pair emission. 
As this is the inverse of the neutrino-pair absorption process, the relevant strength function 5'*^"** can 
be obtained by detailed balance 


S“(A,r) = S“‘*(A)xexp('-TV (24) 

While the absorption strength is independent of temperature by construction, the emission strength 
function depends on temperature by the Boltzmann factor exp {—A/kT) = exp(—(77^, -|- Ep)/kT). 

The ansatz Eq. (23) is a rather simple approximation. For example, it assumes that it is the 
same for all nuclei and it neglects the dependence of the widths on temperature. An alternative 
approach to describe the nuclear deexcitation rate has been suggested by Fuller and Meyer [194]. These 
authors derived analytical expressions for the allowed and forbidden contributions to the emission and 
absorption strength functions at hnite temperature based on the Fermi gas model. The parameters 
in their expressions were adjusted to reproduce results obtained in an independent single-particle shell 
model. We note that the strength functions proposed in Ref. [194] depend on temperature as well as 
on charge and mass number of the nucleus. 

The absorption and emission strengths of the two models are plotted in Fig. 19. The absorption 
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strength functions of Ref. [194] show two distinct peaks for allowed and forbidden transitions, while 
the Gaussian ansatz (Eq. 23) distributes the strength over a wider energy range. While the Gaussian 
ansatz is independent of charge and mass number of the nuclei, the absorption strength of Ref. [194] 
is proportinal to the number of nucleons and increases with temperature as the NSE composition gets 
shifted to heavier nuclei with growing temperature. For example, the average mass number is (A) = 124 
at T = 1.5 MeV using the Lattimer-Swesty Equation of State [71] and (A) = 205 at T = 3 MeV. In 
summary, the two absorption strengths are quite distinct and can be understood as extreme cases 
allowing to explore how sensitive supernova simulations are to an accurate description of the strength. 
Nevertheless it might be desirable that the absorption strength function will be calculated within the 
TQRPA approach for several nuclei and at different temperatures. 

The emission strength function is obtained by multiplying the absorption strength with the appro¬ 
priate Boltzmann factor, see Eq. 24. The differences between the two absorption strength functions lead 
to strong deviations in the emission strength functions. We note that there is also emission strength at 
energies above the thermal average energy due to the thermal population of states at higher energies. 
In the emission decay rate the high-energy tails of the emission strength function get enhanced by an 

factor in the phase space factor, while this factor relatively decreases the contributions from lower 
energies. 

Recently the downwards strength function has been calculated for hxed exciation energies for several 
nuclei within the diagonalization shell model [196]. To translate the downward strength into a thermal 
emission strength, the excitation energy has been adequated with the average thermal energy of a Fermi 
gas model. Hence by construction this approach does not allow for the production of neutrino pairs with 
a sum energy larger than the thermal average energy. Due to the phase space factor such high-energy 
pairs do noticeably contribute to the deexcitation rate. 

The impact of nuclear deexcitation on the supernova evolution has been studied in core-collapse 
supernova simulations of an 11.2 Mq progenitor star [195]. Three simulations have been performed 
which differ only in the treatment of the neutrino pair de-excitation process. In the hrst two simulations 
this process has been implemented following the prescription of Ref. [194] and by adopting the Gaussian 
approximation model to the absorption strength function, as described above. The third simulation 
served as a control study, in which the neutrino-pair nucleear deexcitation has been switched off. 

The most important result of these simulations is that the global quantities such as temperature and 
electron fraction profiles are the same in all simulations. This implies that the neutrino-pair nuclear 
deexcitation process has no impact on the global supernova evolution [195]. This is confirmed in Fig. 
20 which shows that the luminosities of electron neutrinos, arising mainly from electron capture, are 
about 4 orders of magnitude larger than those of heavy-flavor neutrinos during the collapse phase. 

Fig. 20 also shows that the evolution of the luminosities and average energies are the same in 
all three simulations: the luminosity increases from about 10®^ erg s“^ to 10^^ erg s“^ at Ue trapping 
shortly before core bounce caused by the increasing density which enlarges the electron Fermi energy 
and the electron capture rates. Relatedly, the average energies of Ue neutrinos increase from 6 MeV 
to 10 MeV which, however are smaller than those of the neutrinos directly produced by the capture 
process, reflecting the importance of down-scattering by interaction with matter. 

Nuclear deexcitation by neutrino-pair emission produces all neutrino types, while electron capture 
is a pure source of electron neutrinos. Indeed the de-excitation process is the dominating source of Ux 
neutrinos and, to a lesser extent, of electron antineutrinos during collapse (see Fig. 20). Due to the 
charged-current contribution the production of electron neutrino pairs is larger than that of the other 
two flavors. At high densities of order 10^^ g cm“^ also nucleon-nucleon bremsstrahlung becomes a 
significant source of neutrinos other than Uf. which causes the steep rise of the and z/g luminosities 
at times just before bounce in Fig. 20. At densities below 10® g cm“® electron-positron annihilation is 
another relevant neutrino-pair production process. The average energies of the heavy-flavor neutrinos 
increase slightly during the collapse (as temperature increases), but they are noticeably smaller than 
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Figure 20: Evolution of neutrino luminosities and average energies from a core-collapse supernova 
simulation of an 11.2 Mq progenitor star including the production of neutrino pairs from heavy-nuclei 
de-excitations, based on Ref. [194] (green lines) and the Gauss ansatz (see text, red lines), in comparison 
to a simulation that uses identical input physics, but neglects the nuclear deexcitation process (blue 
lines) (from [195]). 
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the values of electron neutrinos produced by electron capture. For densities p > 5 x 10^^ g cm“^ 
neutrino-matter interactions also become relevant for heavy-flavor neutrinos. 

At core bounce where normal matter density is reached, heavy nuclei dissociate into a state of 
homogeneous matter of nucleons and hence the production of neutrino pairs from nuclear deexcitation 
disappears. At shock formation and during the initial shock propagation out of the stellar core, the 
infalling heavy nuclei that hit the expanding shock wave also dissociate. Consequently at the conditions 
behind the expanding shock front, weak processes are determined by interactions with free neutrons 
and protons. Hence, the inclusion of heavy-nuclei de-excitations has no impact on the supernova 
dynamics or the neutrino signal after core bounce, e.g., in terms of the energy loss (as suggested in 
ref. [194]). Although a small fraction of heavy nuclei exist ahead of the expanding bounce shock before 
being dissociated (see the discussion in the preceding subsection), the conditions are such that other 
pair processes dominate over the pair production from nuclear de-excitation. Moreover, the supernova 
dynamics is dominated by charged-current processes on free nucleons behind the bounce shock in the 
dissociated regime. Consequently, neutrino pair heavy-nuclei de-excitations has no impact in the entire 
post-bounce period and the evolution of the neutrino luminosities and average energies in our three 
simulation become identical. 

By studying the radial prohles of the collapsing core at different times during collapse. Ref. [195] 
gives a detailed account of the evolution of the luminosities and average energies of the various neutrino 
flavors and the competition of the different neutrino sources. It turns out that the two approaches 
which have been adopted to describe nuclear deexcitation give qualitatively similar results, they differ, 
however, in details. This is relevant if one is interested in an accurate description of the heavy-flavor 
and electron antineutrinos during collapse. An improved calculation of the absorption strength function 
within a formalism like the Thermal QRPA of Ref. [134] might be useful to this end. 


5 Neutrino-nucleus reactions and supernova nucleosynthesis 

As discussed above, the continuous emission of neutrinos from the protoneutron star drives a low-mass 
outflow (the neutrino-driven wind). Due to the high temperature involved the matter is ejected as free 
protons and neutrons. Upon reaching cooler regions further away from the neutron star, the matter can 
assemble into nuclei. The outcome of this nucleosynthesis depends crucially on the proton-to-neutron 
ratio and can either support the i^p process (if > 0.5) or an r process (Y). < 0.5). We will discuss 
both processes in the two following subsections. Furthermore neutrinos can also contribute to the 
astrophysical production of certain nuclides, in what is called the u process, when they pass through 
the outer shell of the star. The z/ process is discussed in the third subsection. 

The proton-to-neutron ratio is mainly determined by the neutrino and antineutrino charged-current 
reactions with the free neutrons and protons [197]. Assuming that the neutrino luminosities are about 
equal for electron neutrinos and antineutrinos a detailed analysis of these reactions shows that the matter 
is expected to be proton rich if 4A > {{EpJ — where A = 1.29 MeV is the proton-neutron mass 

difference and {EpJ, {E^^) are the average energies of electron antineutrinos and neutrinos, respectively 
[198]. It is important that in the kinematics of the neutrino-nucleon reactions corrections due to the 
nuclear interaction are included. The dominant contribution is due to mean-held corrections, induced by 
the astrophysical environment and accounted for in the Equations of State [199, 200]. Introducing the 
mean-held potentials for protons. Up, and neutrons, Un, the energy of neutrinos is reduced by Un — Up, 
while the energy of antineutrinos is enlargened by the same amount. We note that these corrections 
act like the Q value in electron capture reactions on nuclei which grows with increasing neutron excess 
throttling the capture rate. In fact, for the neutron-rich environment in core-collapse supernovae, e.g. 
shortly after bounce, Un — Up > 0 and as a consequence, the emissivity of neutrinos (antineutrinos) is 
increased (decreased) if the mean-held corrections are taken into account. 
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Figure 21: Comparison of elemental abundances in the mass range Ca to Zn produced in proton-rich en¬ 
vironment (Frohlich et al., [198]) with previous calculations (Thielemann et al., [212]) and observational 
data [213, 214] (adapted from [198]). 


Remembering that the charged-current neutrino reactions on nucleons set the neutron-to-proton 
ratio of the matter ejected in the neutrino-driven wind, the changes of the neutrino and antineutrino 
energies due to the mean-held corrections have an essential impact on the neutrino-driven wind nucle¬ 
osynthesis: the early ejecta are slightly neutron rich, while at later times the ejecta become proton rich 
[ 200 ], 

In Ref. [201] it has been shown that in the outer layers of the protoneutron star a substantial amount 
of light nuclei (deuterons, ^H, ^He, and "^He) is present. The interaction of neutrinos with these light 
nuclei has an impact on the average energies of neutrinos and antineutrinos. The relevant neutrino 
cross sections on the light nuclei have been calculated on the basis of ab-initio few-body methods like 
the hyperspherical harmonics approach [202, 203, 204]. These cross sections are more accurate than 
previous calculations performed on less sophisticated nuclear models (e.g. [9, 205]), but do not differ 
too dramatically. 

5.1 The up process 

Supernova simulations indicate that during some period after bounce the hot matter ejected from the 
surface of the freshly born neutron star is proton-rich (i.e. Fj. > 0.5) (e.g. [50, 206, 207, 208, 8]. The 
alpha rich freeze-out of such proton-rich matter favors the production of a nuclei (mainly ®®Ni) with 
some free protons left [209]. We note that this freeze-out also results in enhanced abundances of selected 
nuclei in the Ca-Fe mass range bringing them into better agreement with observation [198, 210, 211] 
(see Fig. 21). 

Heavier nuclides can be synthesized from the freeze-out abundance distribution by subsequent proton 
captures, competing with /3+ decays. We note that this reaction sequence is realized in the so-called 
rapid-proton capture process (i.e. explosive hydrogen burning on the surface of neutron stars in X-ray 
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Figure 22: (left) Evolution of the abundances of neutrons, protons, alpha-particles, and ^®Ni. The 
solid (dashed) lines display the nucleosynthesis results which include (omit) neutrino and antineutrino 
absorption interactions after nuclei are formed. The abscissa measures the time since the onset of the 
supernova explosion, (right) Isotopic abundances (relative to solar) obtained within the supernova up 
process nucleosynthesis. The hlled circles represent calculations where (anti)neutrino absorption reac¬ 
tions are included in the nucleosynthesis while for the open circles neutrino interactions are neglected. 
The two panels base on a nucleosynthesis trajectory from the model B07 from Ref. [198] (from [210]). 


bursts) [215]. However, in this scenario the mass flow to heavier nuclides is strongly hampered by the 
increasing Coulomb barrier of the produced elements and, in particular, by the so-called waiting point 
nuclei. These are a nuclei like ^®Ni, ®^Ge, ®®Se... which have relatively long (3 half lives and for which 
the proton capture is strongly hindered due to the small or negative proton binding energies of the hnal 
nuclei (e.g. ®^Cu, ®®As, ®®Br, ...). In the case of unbound daughters it is suggested that the mass flow 
might overcome the waiting point nuclei by two-proton captures [216]. 

In contrast to the nucleosynthesis in X-ray bursts, the process occurs in the supernova environment 
in the presence of extremely intense neutrino fluxes which influences and alters the matter flow to 
heavier nuclei substantially. While the energy of supernova z/g neutrinos is too small to induce sizable 
reaction rates on ~ Z nuclei, this is different for antineutrinos that are captured in a typical time of 
a few seconds in the conditions of the hot neutrino bubble, both on protons and nuclei, at the distances 
at which nuclei form (~ 1000 km). This time scale is much shorter than the beta-decay half-life of 
the most abundant heavy nuclei (e.g. ®®Ni, ^'^Ge), which serve as obstacles in the mass flow to heavier 
nuclei if neutrino interactions are not included. As protons are more abundant than heavy nuclei, 
antineutrino capture occurs predominantly on protons, causing a steady supply of free neutrons with 
a residual density of lO^'^-lO^® cm“^ for several seconds, when temperatures are in the range 1-3 GK 
(see left panel in Fig. 22). The neutrons produced via antineutrino absorption on protons can easily 
be captured by neutron-dehcient 77 ~ Z nuclei (for example ®^Ge), which have large neutron capture 
cross sections. The amount of nuclei with A > 64 produced is then directly proportional to the number 
of antineutrinos captured. While proton capture, (p, 7 ), on ®^Ge takes too long, the (n,p) reaction 
dominates (with a lifetime of 0.25 s at a temperature of 2 GK), permitting the matter flow to continue 
to nuclei heavier than ®^Ge via subsequent proton capture up to the mass range A ~ 80 — 100 (for 
a mass-flow diagram see, for example. Ref. [217]). This nucleosynthesis process, operating in proton- 
rich supernova environment, is called the z^p-process [210]. The abundance distribution obtained by 
z/p-process nucleosynthesis is compared in the right panel of Fig. 22 to a simulation in which the effects 
of (anti)neutrino captures are switched off. Huther has explored the influence of the various neutrino- 
induced neutral- and charged-current reactions on the up process abundances. He hnds that the impact 
of the neutrino and antineutrino reactions on nuclei is negligible [218]. 
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How far the massflow within the up process can proceed strongly depends on the environment 
conditions, most noticeable on the valne of the matter. Obvionsly the larger Yg, the larger the 
abnndance of free protons which can be transformed into nentrons by antinentrino absorption, and the 
faster the bridging of the waiting point nnclei like ®^Ge by {n,p) reactions. Other important parameters 
are the location (radius) of the matter during the formation of nuclei and the ejection velocity which 
both influence the antineutrino fluence which the ejected matter experiences. A location closer to the 
surface of the proto-neutron star and/or a slow ejection velocity leads to an extended antineutrino 
exposure that allows for an increased production of heavy elements [219, 220]. The sensitivity of the 
z/p-process on parameters like the T/ value has been explored in [210, 217, 221] These studies show that 
nuclei heavier than A = 64 are only produced for Yj. > 0.5, showing a very strong dependence on Tj. 
in the range 0.5-0. 6 . A clear increase in the production of the light p-nuclei, and ®®’®®Ru, is 

observed as Tj. gets larger. Thus the up process offers an explanation for the production of these light 
p-nuclei, which was unknown before. 

The late-time abundance of free neutrons depends sensitively on the Ue luminosities and spectra. As 
has been shown recently [94, 93, 95], the Ue spectrum is affected by collective neutrino flavor oscillations 
which are expected to occur in the high-neutrino-density environment surrounding the neutron star [93] 
and can swap the Ue and u^^r spectra above a certain split energy of order 15-20 MeV, where the split 
energy depends on the relative fluxes of u^ and u^ and on the neutrino mass hirarchy [ 222 ]. Assuming 
such a swap scenario for the case of the normal neutrino mass hierarchy (see left part of Fig. 23), Ref. 
[223] has studied the impact of collective neutrino flavor oscillations on the up process nucleosynthesis 
adopting anti-neutrino spectra from recent supernova simulations [90]. As in general, Ux neutrinos 
have larger mean energies and a larger high-energy tail than electron antineutrinos, the spectrum swap 
induced by collective neutrino oscillations increases the neutron production rate by antineutrino capture 
on protons during the up process. Ref. [223] estimates this increase to vary by a factor between 1.5 to 
1.25 depending on a split energy between 0 and 25 MeV. Obviously the enhanced production of free 
neutrons boosts the matter flow to heavier nuclides and results in larger up process abundances. This 
is demonstrated in the right part of Fig. 23 which compares the z/p-process abundances calculated with 
a neutron production rate increased by a factor of 1.4 due to collective neutrino oscillations and the 
standard case without such oscillations. The production of the light p-nuclides and ®®’®®Ru is 

enhanced by a factor 2-3. The largest production increase is observed for nuclides A > 96 which is 
traced back to a faster (n,p) reaction time-scale for the N = 50 nucleus ®®Pd which acts like a ’seed’ 
for the production of nuclides with A > 96 [224]. The faster (n,p) reactions also have an effect on 
the abundances of waiting point nuclei like ®‘^Ge or ®®Se. These nuclei are less of an obstacle in the 
mass flow in the calculations which consider collective neutrino oscillations and hence less matter is 
accumulated at the waiting points. More recently, Wu et al [97] have performed a detailed analysis of 
the impact of collective oscillations on neutrino-wind nucleosynthesis accounting for the full temporal 
variation of neutrino spectra and matter density profiles. 

The up process constitutes a mass flow from seed nuclei, mainly ®®Ni, to heavier nuclides by (p, 7 ) 
reactions, supplemented by (n, p) or (3^ reactions running through many nuclei at the proton-deficient 
side of the nuclear chart. Many of the (3^ half-lives, for the ground states, are experimentally known. 
However, their potential modihcations for the astrophysical environment with hnite temperature as well 
as the rates for the proton capture and (rz,, p) reactions have to be theoretically modelled, which is usually 
performed on the basis of the statistical model. A particularly important ingredient in such statistical 
model calculations are the nuclear masses, which dehne the reaction Q values and thus the competition 
between reactions and their inverse. Nuclear mass measurements have benefited tremendously in recent 
years by the establishment of nuclear Penning traps or by storage ring experiments at radioactive ion- 
beam facilities which also allowed the determination of masses of z/p-process nuclei with unprecedented 
precision. With relevance to the up process such measurements have been performed for nuclei in the 
mass range A ~ 90 at the SHIPTRAP at GSI and JYFLTRAP in Jyvaskylla [225, 226], for nuclei 
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Figure 23: Left: (upper panel) Pg and Vx spectra /p^ and /p^ based on the simulations of ref. [90] as 
given in table 1 of ref. [217]. The Pg spectra have been arbitrarily normalized to one while the Vx spectra 
are normalized to 1.3 to keep the relative ratio of the number luminosities, (lower panel): modified 
spectra /p^ and /p^, including the effects of collective neutrino oscillations, which induce a swap of the 
Pg and Vx spectra for energies above Eg. Right: Abundances of vp nucleosynthesis calculations with 
(open circles) and without (full circles) considering the effects of collective neutrino flavor oscillations 
(normalized to solar). The lower panel shows the ratios for the two nucleosynthesis studies as function 
of mass number (from [223]). 
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around A = 64 at the CSRe storage ring in Lanzhou [227], for proton-dehcient nuclei in the mass range 
A ~ 90 at the Canadian Penning trap [228]. The impact of these mass measurements on the up process 
has been explored in Refs. [224, 225]. Studies which investigated the sensitivity of the up process 
abundances on the uncertainties in the {n,p) rates have been reported in [219, 229]. The results are 
similar to those found for the effect of the collective neutrino oscillations: faster (n,p) rates spead up 
the mass flow to heavy nuclides and hence increase the abundances for nuclides with A > 96, while 
they decrease the abundances for nuclides in the mass range A ~ 64 — 96. 

5.2 The r-process 

More than 50 years ago the astrophysical r-process (short for rapid neutron capture process) was 
introduced to explain the origin of the transactinides and of half of the nuclides heavier than iron 
[19, 230]. The process is characterized by fast neutron captures, in comparison to the competing (3 
decays, which, in turn, requires astrophysical objects with extreme neutron densities for the r process to 
operate [108]. On the nuclear chart the r-process path runs through nuclei with extreme neutron excess 
far off the valley of stability. Most of these nuclei have yet not been produced in the laboratory and their 
properties, which are required input to r-process simulations, have to be modelled theoretically. The 
most relevant nuclear input are masses (reaction Q values) and f3 half lives which influence the r-process 
path and the duration of the process, respectively [108]. Particularly important for the r-process mass 
flow to heavier nuclei are the waiting point nuclei associated with the magic neutron numbers N = 82 
and 126 (and possibly 184), which are also related to the peaks observed in the r-process abundance 
distributions around mass numbers A ~ 130 and 195. 

Decisive progress towards experimental information for r-process nuclei is expected from current and 
future radioactive ion-beam facilities like RIBF at RIKEN, FAIR at GSI and FRIB at NSCL/MSU. 
From these facilities we expect not only the direct determination of properties of r-process nuclei, but 
also stringent constraints to improve the theoretical description of such neutron-rich nuclei which are out 
of experimental reach even at these advanced facilities. We have witnessed such indirect progress due 
to data for neutron-rich nuclei for both, the global description of masses and the prediction of half-lives. 
The advances in experimental mass measurements are reviewed in Ref. [231]; recent mass evaluations 
can be found in [232, 233] (see also [234]). These data help to improve the predictions of nuclear mass 
models like the microscopic-macroscopic approach for many years pioneered by Moller [235, 236] (see 
also [237]), the Extended Thomas Fermi Model with Strutinski Integral [238], the shell-model guided 
approach by Duflo and Zuker [239] or studies based on the Hartree-Fock-Bogoliubov model [240] . Half- 
life measurements of r-process nuclei at or in the vicinity of the magic neutron numbers 77 = 50 and 
82 have been reported in Refs. [241, 242, 243, 244] and [245, 246, 247], respectively. First half-life 
measurements towards r-process nuclei at the N = 126 shell closure have been performed at GSI [248], 
which are useful constraints for theoretical half-life predictions in this until yet unexplored part of the 
nuclear chart (e.g. [249, 250, 251] and the calculations of Borzov as quoted in [248]). 

One of the fundamental questions in astrophysics is: where does the r-process operate? Although no 
conclusive answer can yet be given, crucial clues come from abundance observations in so-called metal- 
poor stars, i.e. stars with Fe/H ratios which are signihcantly (often by orders of magnitude) smaller 
than the solar Fe/H ratio. While hydrogen has mainly been produced by Big Bang nucleosynthesis, iron 
has gradually been synthesized by both types of supernovae (core-collapse and thermonuclear) during 
the history of the galaxy. The Fe/H abundance ratio can thus serve as a proxy of age and stars with 
Fe/H ratios much smaller than solar are supposed to be very old stars. Due to their abundance pattern, 
metal-poor stars can be classihed roughly into two categories: So-called r-H stars [252] exhibit relative 
abundance patterns which agree nearly perfectly with the solar r-process abundances for elements with 
charge numbers Z > 50 [253, 254, 255]. Some other stars [256, 257] are depleted in elements with 
Z > 50, but enriched in lighter elements like Sr, Y and Zr, with the star HD 122563 being a famous 
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Figure 24: R-process abundances calculated for a neutron-star merger trajectory and different mass 
models (from [258]). 


example [256]. One often denotes the process whose imprint is observed in stars like HD 122563 the 
’weak r-process’, while the ’main r-process’ is related to the robust abundance pattern observed in the 
r-II stars. 

It is tempting to identify potential astrophysical sites for the r-process. The main r-process might 
be attributed to neutron star mergers. Recent simulations [259, 260] show that the extreme neutron 
richness of the matter dynamically ejected during the merger {Yg ~ 0.05) leads to fast mass flow to 
nuclei up to mass range A ~ 280 and the occurence of fission cycling. In fact, the second r-process peak 
at H ~ 130 is produced by hssion fragments in this scenario. Hence fission rates and yield distributions, 
including the production of free neutrons during the hssion process and by /3-delayed neutron emission 
of hssion yields, are crucial ingredients in the simulations [261, 262]. Fig. 24 shows the abundance 
distribution obtained from r-process simulations based on neutron-star merger trajectories of Ref. [260] 
and using the hssion fragment yields calculated on the basis of the ABLA code [263] which has been 
carefully adjusted against experimental hssion data [111]. In fact the calculated abundance distributions 
reproduce quite nicely the solar r-process pattern in a quite robust way, i.e. basically independent on 
the nuclear mass models and the details of the merger trajectories [258]. Despite these intriguing 
results it is still an open question whether neutron-star mergers have already operated with sufficient 
frequencies at the early time of the galaxy to explain the observed r-process abundances in the oldest 
stars [264, 265, 266, 267, 268]. Furthermore, the impact of neutrinos in relativistic merger models 
requires further investigation [269]. Alternative sites for the r-process which are expected to operate 
mainly at low metallicities are jets from magnetorotational supernova [270] or an r-process in the He- 
shell of core-collapse supernovae [271]. In the latter scenario, the required neutrons are produced by 
charged-current reactions on ^He via "^He(t'e, where the relevant cross sections are given in [272] 

based on an ab-initio calculation within the hyperspherical harmonics approach [273]. 

The weak r-process has been suggested to operate in the neutrino-driven wind developing in core¬ 
collapse supernovae after bounce. This scenario has actually been the favorite site for the r-process for 
many years [275, 276]. To achieve the neutron-to-seed ratio n/s necessary to produce the third r-process 
peak at A ~ 195 and the transactinides after freeze-out of charged-particle reactions, n/s > 150, puts 
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Figure 25: Mass-integrated abundances for the neutrino-driven wind of a core-collapse supernova 
(adapted from [274]). 


stringent constraints on the wind properties. Such high ratios require either large entropies, small 17 
value or fast velocities of the ejected matter [277, 278]. Recent supernova simulations imply that such 
conditions are not met in the neutrino-driven scenario making the operation of the ’main r-process’ 
unlikely in this scenario. However, the conditions might be sufficient to produce the light r-process 
elements up to Sr, Y, and Zr [279, 274]. This is demonstrated in Fig. 25 which shows the mass- 
integrated abundance distribution (upper panel) and elemental abundances (lower panel) obtained from 
a core-collapse supernova simulation of a 11.2Mq progenitor star [274] which followed the post-bounce 
evolution up to 9 s after bounce. The calculation adopted an Equation of State which consistently 
provided the mean-held corrections needed for the calculation of the charged-current neutrino reactions 
on nucleons which determine the Ye value of the ejected matter (as discussed above). The ejected 
matter was found to be neutron rich at early times (Yg ~ 0.48 for the hrst 3 seconds after bounce). At 
later times the matter becomes increasingly enriched in protons reaching Yg = 0.58 after 9s. The early 
neutron-rich ejecta produce isotopes of the elements Sr to Mo (Z = 38 — 42), while nuclei with Z > 42 
(or A > 92) are produced by the i/p process in the later proton-rich ejecta. However, their production 
is quite inefficient due to low antineutrino luminosities (equivalently neutron production rates) at late 
times. Ref. [274] also hnds relatively small amounts of matter ejected at late times. Thus, the z/p 
process does not contribute much to the mass-integrated abundances. It is very interesting to note 
that the elemental abundances for the elements Sr, Y, Zr, Nb, Mo reproduce the abundance pattern 
observed in the star HD 122563 quite well (see lower panel in Fig. 25). However, the calculation clearly 
underproduces the observed abundances for the elements with Z > 42, including silver and palladium. 
The production site of these elements is yet an open question [274]. Active-sterile neutrino oscillations, 
that assume a sterile neutrino mass ~ 1 eV, help to achieve conditions suitable for the productions of 
elements in the range Z = 42-50 [280] 

There have been several suggestions how neutrino-nucleus reactions might alter the r-process abun¬ 
dance in the neutrino-driven wind scenario (for a review see [281]). However, the impact of the re¬ 
actions has mainly been explored chosing parametrizations of the wind which made it the site of 
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the main r-process. The studies include r-process calculations including neutrino-nucleus reactions 
[282, 283, 284, 285, 286], the exploration of the role of neutrino-induced hssion [106, 107, 287] and 
potential hngerprints left by neutrino-induced spallation on the final r-process abundances [288, 289]. 
The recent calculations by Huther [218] indicate that neutrino-nucleus reactions have no noticeable 
impact on the abundance distribution of the weak r-process [274], A potential exception is one of the 
variants of the so-called ’a effect’, proposed in Ref. [290]. In very strong neutrino fluxes neutral-current 
reactions on “^He can trigger two-body reaction sequences (e.g.^He(z/, z/'p)^H(^He, 7 )^Li(^He, 7 )^^B and 
‘^He(i/, z/'? 7 ,)^He('^He, 7 )^iBe('^He, 7 )^^C) which compete with the three-body reaction a + a + n ^ ®Be 
which has been initially considered to bridge the mass gaps at A = 5 and 8 in the neutrino-driven wind 
r-process [275]. As a consequence of the more effective bridging the mass flow to heavier nuclei is faster, 
reducing in turn the neutron-to-seed ratio and thus throttling the subsequent r-process. 

5.3 The TJ process 

When neutrinos, produced in the hot supernova core, pass through the outer shells of the star, they can 
induce nuclear reactions and in this way contribute to the synthesis of elements (the i/-process) [291, 292]. 
As pointed out by Woosley et al. [9], the nuclides ^^B and are produced by (i/, z/'u) and 
reactions on the quite abundant nuclei and ^‘^Ne. These reactions are dominantly induced by and 
Vr neutrinos and their antiparticles due to their larger average energies than z/g and z7g neutrinos. Also 
the synthesis of the odd-odd nuclides ^^®La and ^^‘^Ta are attributed to the v process. However, they 
are produced mainly by the charged-current reaction ^^®Ba(z/g, e“)^^®La and ^®°Hf(z/g, e“)^®°Ta [293]. 
Hence, the z/-process is potentially sensitive to the spectra and luminosity of z/g and Vx neutrinos, which 
are the neutrino types not observed from SN1987A. 

Neutrino nucleosynthesis studies are quite evolved requiring state-of-the-art stellar models with an 
extensive nuclear network [186, 294]. In the first step stellar evolution and nucleosynthesis is followed 
from the initial hydrogen burning up to the presupernova models. The post-supernova treatment then 
includes the passage of a neutrino flux through the outer layers of the star and of the shock wave 
which heats the material and also induces noticeable nucleosynthesis, mainly by photodissociation (the 
7 process [295]). Modelling of the shock heating is quite essential as the associated 7 process destroys 
many of the daughter nuclides previously produced by neutrino nucleosynthesis. 

The production of ^^®La and ^®°Ta by the v process is caused by charged-current (z/g,e“) reactions 
on ^^®Ba and ^®°Hf, respectivley. At the energies involved the cross sections are dominated by the 
low-energy tail of the GT_ strength distribution. As large-scale shell model calculations are yet not 
feasible for these relatively heavy nuclei, the cross sections had to be calculated within the RPA model 
for both allowed and forbidden transitions [9, 294]. As a major improvement the GT_ strengths on 
^^®Ba and ^®°Hf below the particle thresholds was measured at the RGNP in Osaka using the (^He, t) 
charge-exchange reaction [296]. The left panel of Fig. 26 compares the running sum of the measured 
B(GT_) strength with the RPA prediction. Note that for the z/-process production of ^^^La and ^®°Hf 
one ist only interested in the strength below the neutron threshold (indicated by Sn in Fig. 26) as the 
strength above connects to states which will subsequently decay by neutron emission and hence not 
contribute to the ^^®La and ^®°Hf production. One hnds that the RPA prediction are slightly smaller 
than the B(GT_) data for ^^®Ba up to the neutron threshold, but underestimates the ^^‘^Hf strength 
by nearly a factor of 3. The inversion of the B(GT_) strength into cross section is complicated by the 
fact that for both nuclei the proton decay channel opens below the neutron channel. The respective 
corrections have been derived from branching ratios obtained from the statistical model. Furthermore 
adding forbidden contributions to the cross sections derived from RPA calculations, one finds that the 
^^®Ba(z/e, e“)^^®La and ^®°Hf(z/e, e“)^®‘^Ta, calculated for a supernova spectrum for z/g neutrinos with 
temperature T = 4 MeV and zero chemical potential, are about 25% and 30% larger than estimated 
solely on the basis of allowed transitions. The cross sections are shown, together with other relevant 
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Figure 26: Summed B(GT) strength in and ^®°Ta as a function of excitation energy up to the 

respective neutron emission thresholds. Solid lines: experimental data; dashed lines: RPA calculation 
[294] (from [296]). 
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Figure 27: (Left) Electron neutrino charged-current cross sections on for cascades up to two 

particle emission, (right) Production of ^^®La in a 25 Mq star. The hgure shows the production of ^^®La 
without neutrinos (dotted line), with the charged-current reaction (dashed), with the 

neutral current reaction ^^^La.{u,u'n) (dashed-dotted) and with both reactions (solid). Additionally, 
the supernova shock temperature, the pre-supernova density, and the neutrino energy fluence 
assuming a total neutrino energy of 3 x 10^^ erg are given (from [294]). 


cross sections for the study of neutrino nucleosynthesis in the left panel of Fig. 27. 

Incorporated into stellar models [186, 294] one hnds indeed that ^^®La is being produced by charged- 
current reactions on ^^®Ba, This is demonstrated in Fig. 27 which shows the ^^®La production in the 
oxygen/neon shell of a 25 Mq star [294]. The key is the enhancement of ^^®Ba by an s-process prior 
to the supernova explosion. A ^^®La abundance is basically not existing prior to the passing of the 
supernova neutrinos and is the mainly produced by the (z/e,e“) reaction on ^^®Ba (right panel). In 
principle, ^^®La can also be made by the neutral-current neutron-spallation reaction (i/, z/'n), but its 
contribution to the total ^^®La abundance is insignihcant (right panel). The 7 process produces also 
some ^^^La when the shock passages through the oxygen/neon layer. But its contribution to the total 
^^®La yield is about an order of magnitude less than the charged-current contribution (right panel). In 
fact, ^^®La appears to be co-produced with nuclides like and ^^Mg in massive stars. Using these 
nuclei as tracers for the contribution of massive stars to the solar abundance [186, 294, 296], the observed 
^^®La abundance appears to be mainly due to the u process. 

A comparison of the ^^‘^Ta abundance calculated in stellar models is complicated by the fact that, on 
the time scales of the u process, ^®°Ta can exists in two states: the J = 9~ isomeric state at excitation 
energy = 75.3 keV with a half life of more than 10^® years and the ground state which decays 
with a half live of 8.15 hours. Also the B(GT_) data cannot distinguish between the contributions 
to these two states. Furthermore in the hnite temperature astrophysical environment the two states 
couple via excitation of states at intermediate energies [297, 298]. To treat this coupling, Hayakawa et 
al. have recently proposed a model [299, 300] in which they distinguish between states build on the 
ground and isomeric states (guided by data of Refs. [301, 302]) and assume these two band structures 
couple only weakly to each other. They then followed the decoupling of the two-band structure from 
thermal equilibrium to freeze-out in a time-dependent approach assuming an exponential decrease of 
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Figure 28: (left) Nuclear reaction flow path for the production of light elements induced by neutrino 
reactions on during a supernova explosion, (right) Dependence of the abundance ratio ^Li/^^B on 
the neutrino mixing angle ^13 for normal and inverted neutrino mass hierarchies. The calculations have 
been performed for two sets of neutrino cross sections and branching ratios (from [128]). 


temperature with time. In their calculation 39% of the ^®°Ta produced in the v process survives in 
the long-lived isomeric state. A similar branching ratio has been estimated earlier by Mohr al. [303]. 
Combining these branching ratios with the total ^®°Ta u process abundance as given in [296] makes the 
u process a potential ^®°Ta production site. We note, however, that also s-process nucleosynthesis has 
been suggested to produce nearly 100% of the solar ^®°Ta abundance [304]. 

Following the work of Ref. [9], several studies have conhrmed that the u process can produce the 
nuclides and ^®F in solar abundances by (z/, v'n) and (z^, v'p) reactions on and ^°Ne. These 
studies considered improved stellar simulations [186, 305, 294] and more elaborate neutrino-induced 
spallation cross sections (e.g. [294, 306, 307]). The shell model has been shown to reproduce the GT 
distributions for nuclei in the p shell quite well [306]. Some ^Li is produced by neutrino-induced reactions 
on ^He in the He layer of the star (the (z/, z/'pa) reaction on can in principle also produce ^Li in the 
oxygen/carbon layers, but the respective branching ratio is too small to be significant [18]). 

Very recently Sieverding et al. have performed studies of neutrino-induced nucleosynthesis in su¬ 
pernovae for stars with solar metallicity and initial main sequence masses between 15 and 40 Mq [308] . 
These studies improve the previous investigations i) by using a global set of partial differential cross 
sections for neutrino-induced charged- and neutral-current reactions on nuclei with charge numbers 
Z < 78 and ii) by considering modern supernova neutrino spectra which are noticeably shifted to lower 
energies (i.e. = 2.8 MeV, = 4.0 MeV assuming a Fermi-Dirac neutrino distribution with 

zero chemical potential) compared to spectra previously adopted and discussed above (e.g. [9, 294]). 
Mainly due to the change in neutrino spectra, Sieverding finds slightly smaller abundances for ^Li, ^^B, 
^^®La and ^®°Ta, however, he conhrms the production of these nuclides by neutrino nucleosynthesis [308] 
(see Table 4). 

Interestingly neutrino-induced reactions noticeably contribute to the production of ^®F in more 
massive stars with M > 20 Mq, while they are negligible in lower-mass stars [308]. This is related to 
the fact that ^®F is produced at two different stellar sites: in the C-0 layer by the reaction sequence 
^®0(p,a)^®N(a, 7 ), which is not sensitive to neutrino processes, and in the 0-Ne layer by neutrino- 
induced spallation from the abundant ^°Ne. As in lower-mass stars the 0-Ne layer is less massive 
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Table 4: Production factors relative to solar abundances normalized to production. The table 
shows results from calculations without neutrinos and with neutrinos assuming recent neutrino energy 
spectra. The highest temperature values correspond to the neutrino spectra adopted in [294]. 


Nucleus 

no u 

T 

I'e 

T- 

T 

= 2.8 MeV 
= 4.0 MeV 
= 4.0 MeV 

T 

I'e 

T- 

T 

= 4.0 MeV 
= 4.0 MeV 
= 6.0 MeV 

15 Mq 

star 





^Li 

0.001 


0.28 


2.54 

i^B 

0.007 


1.43 


6.13 

19p 

1.02 


1.14 


1.31 

i38La 

0.07 


0.67 


1.18 

i«0Ta 

0.07 


1.14 


1.81 

25 Mq 

star 





^Li 

0.0005 


0.11 


0.55 

i^B 

0.003 


0.80 


2.61 

19p 

0.06 


0.24 


0.43 

i38La 

0.03 


0.63 


1.14 

i8°Ta 

0.14 


1.80 


2.81 


(about 1 Mq in the 15 Mq star) than in more massive stars (about 3 Mq in the 25 Mq star) the 
relative importance of neutrino-induced reactions for the total abundance depends strongly on the 
progenitor mass and the stellar evolution. 

The studies of Ref. [308] hnd that neutrino-induced reactions, either directly or indirectly by provid¬ 
ing an enhanced abundance of light particles, noticeably contribute to the production of the radioactive 
nuclides ^^Na and ^®A1, which are both prime candidates for gamma-ray astronomy. However, the stud¬ 
ies do not hnd signihcant production of two other candidates, "^^Ti and ®°Fe, due to neutrino-induced 
reactions. 

The neutrino-induced reactions on ^He (in the He layer of the star) and on (in the oxygen/carbon 
layer) initiate a small nuclear network, which produces a few light elements, in particular ^Li and 
(see left panel of Fig. 28). The main path for the ^Li production is by ^He(i/, z/'p)^H(q!, 7 )^Li and 
'^He(i/, z/'u)^He(a, 7 )^Be(e“, r'e)^Li, while the ^^B production is mainly due to the neutral-current (z/, u'n) 
and (z/, z/'p) reactions on However, in both cases also charged-current reactions (^He(z/e, e“p)^He 
and ^He(z/e, e+n)^H in the case of "^He, and ^^C(z/e, e“p)^^B and ^^C(z/e, e+n)^^C for ^^C) contribute to 
the production of ^Li and ^^B, respectively. This makes the abundance of these two nuclides sensitive 
to neutrino oscillations. Due to the expected hierarchy of average neutrino energies < (E^^), 

neutrino oscillations are expected to increase the average z/g energy and consequently also the charged- 
current cross section induced by supernova neutrinos, as pointed out by Kajino and collaborators, this 
makes the ratio of ^Li and ^^B sensitive to the 613 mixing angle and to the mass hierarchy [307, 310, 311] 
as is demnostrated in the right panel of Fig. 28. Despite this intriguing sensitivity, an accurate derivation 
of the ^Li/^^B abundance ratio requires reliable stellar model calculations and neutrino and nuclear 
cross sections (see [312], but must also consider the production of the elements from other astrophysical 
sources; ^Li is, for example, also produced by Big Bang nucleosynthesis [313]. 

Accounting for neutrino oscillations requires to consider the effect of collective neutrino oscillations. 
This has been recently done in ref. [97] based on neutrino spectra obtained in Boltzmann transport 
simulations [314]. Their results show that collective neutrino oscillations increase the average energy 
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of Vf. neutrinos assuming inverted hierarchy. Hence, the production of ^^^La and ^®°Ta, which are 
both generated by charged-current (z/e,e“) reactions, is enhanced. ^Li and are produced in regions 
beyond the MSW resonance. Here the situation is different for normal or inverted hierarchy. In the 
case of normal hierarchy the increased average energy of v^. neutrinos results in an enhancement of the 
production of ^Li by the reaction ^He(z/e, e“p)^He(Q;, 7 )^Be(e’''z/e)^Li. For the case of inverted hierarchy 
the production of ^^B is enhanced by the reaction chain '^He(Pe, e'''n)^H(a, 7 )^Li(a, 7 )^^B. 

Woosley et al. also suggested that some medium-mass nuclei, most notably and ®®Mn, might 
be produced by neutral-current {u, u'p) reactions on ®^Fe and ^®Ni with subsequent beta decays of the 
daughters ®^Mn and ^®Co, respectively [9]. A recent calculation of the nucleosynthesis in a 25 Mq star, 
based on the stellar evolution of the star performed by Heger and for the first time taking the full set of 
neutrino-induced partial cross sections ([218], see above) into account, indeed confirms this suggestion 
showing that a sizable part of the abundances for V, ^^Mn and ^®Co can be attributed to the f process 
[315]. However, the production of these elements occurs close to the supernova core and therefore the 
amount, which might be ejected, will depend on the matter fallback during the onset of the supernova 
explosion requiring a more detailed study. 


6 Detecting supernova neutrinos 

The observation of neutrinos from supernova SN1987A by the earthbound detectors Kamiokande [1] and 
IMB [316] has conhrmed and advanced the understanding of core-collapse supernovae. A similar boost 
is expected from the observation of the next near-by supernova which is likely to test the predictions 
of supernova models concerning the neutrino spectra for the different flavors, including the noticeable 
neutrinoburst signal in electron neutrinos, originating from electron capture on protons set free by the 
shock, and the expected hierarchy in the average neutrino energies between the different flavors {{E^J < 
{EpJ < {Eu^}) reflecting the differences in the interaction of the neutrino flavors with the surrounding 
matter in the supernova core. To reach this goal, several supernova detectors are operational, with 
some more being in preparation or proposed. 

Obviously translating the event rates of supernova neutrinos observed in the detectors requires a 
detailed understanding of these detectors, including the knowledge of the cross sections for the neutrino 
interaction with the detector materials. For a recent overview on the detectors and their scheme for 
observing supernova neutrinos the reader is refered to Ref. [191]. In accordance with the theme of the 
present review we will focus in the following on recent advances in describing neutrino-induced cross 
sections for those nuclei which serve as detector materials in current and future supernova neutrino 
detectors. Much of the early work on this topic has been described in [18]. We will briefly summarize 
these studies and discuss in some details the more recent calculations. 

Besides neutrino-induced inelastic reactions on nuclei - as we will discuss in the following - supernova 
neutrinos can also be detected by their interaction with electrons and by elastic scattering on nuclei 
in the detector. The cross sections for the interaction of supernova neutrinos with electrons is well 
understood [317], but small. However, it has the advantage that the direction of the neutrino can 
be reconstructed if the electron track after the scatter can be reconstructed (like in Water Cerenkov 
detectors or liquid argon time-projection chambers) [191]. Elastic neutrino scattering on nuclei must 
be detected by the nuclear recoil. This can be a viable scheme for protons [318], but for nuclei the 
corresponding recoil energies are too small [319, 320, 191]. The cross section for this coherent process 
is well known and has been discussed above, see eq. (1). 

Table 2 in Ref. [191] lists the present and future supernova neutrino detectors. Their main material 
are liquid scintillator (C„H 2 n), water, lead or liquid argon. Hence we will in the following discuss in 
turn reactions on ^°Ar and ^°®Pb. Finally we add some remarks on neutrino-induced reactions 

for selected molybdenium and cadmium isotopes motivated by the ability of the MOON and COBRA 
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detectors to potentially also observe supernova neutrinos. 


6.1 Carbon 

The nucleus is an essential part of liquid scintillator detectors like Borexino, KamLAND, Mini- 
BooNe, LVD, Baksan. Electron neutrinos and antineutrinos can interact with by charged-current 
reactions (i.e. e“)^^N and ^^C(z/e, e’'')^^B), but for the relatively low-energy supernova neutrinos 

with {E^^) ^ 10 MeV and (Ep^) ^ 15 MeV both reactions are hindered due to the large Q values 
with 17.34 MeV for Ue and 14.39 MeV for Ug. The reactions leading to particle-bound states in the 
hnal nucleus can be tagged by the beta decays of the and ^^B ground states. The matrix ele¬ 
ment governing the neutrino-induced excitations of either ground state can be readily derived from 
the experimental value of the beta decay. We note that in RPA calculations of the and 

^^C(z/e, e’'')^^B inclusive cross sections the ground state transition is adjusted to the beta-decay data. 
Then the calculated e“)^^N cross section for electron neutrinos, produced by muon decay at rest 

and with somewhat higher energies than supernova z/g neutrinos, agree quite well with the data from 
the KARMEN and LSND collaborations (see [18, 321]). 

Recently shell model calculations have been performed for charged- and neutral-current neutrino 
reactions on using two different residual interactions in a (0-|-2)ha; model space [306]. Two results 
of these studies are quite interesting. At hrst, the B(GT) strength is noticeably reduced if the model 
space is extended from the pure p shell (Ohw) to the ( 0 - 1 - 2 )hen requiring for the latter case nearly no 
quenching of the GT strength anymore. Furthermore, the various neutrino-induced cross sections agree 
quite well with those obtained by RPA calculations (e.g. [322]). Shell model calculations performed in 
model spaces larger than the p shell, but employing different residual interactions than Ref. [306] , have 
also been reported in [126, 127] again showing the importance of correlations beyond the RPA approach 
to describe the neutrino-induced transitions to the T = 1 triad in the A = 12 nuclei. 

The detectors can observe neutral-current reactions, triggered by all neutrino flavors, by detecting 
the 7 ray following the {u, v') excitation of the isospin T = 1 state at excitation energy E^ = 15.11 MeV. 
This state is the Isobaric Analog State of the ^^B and ground states. The Gamow-Teller matrix 
element needed to calculate the ^^G(z/, z/')^^G(15.11) cross section can be derived from the beta-decay 
of the ^^B and ground states after an appropriate isospin rotation or directly from the measured 
Ml strength between these states determined in inelastic electron scattering. 

Due to the high thresholds for all neutrino reactions on ^^G the major detection channel is Ue 
absorption on free protons for liquid scintillator detectors. The respective event rate is estimated to be 
about an order of magnitude larger than the charged-current reactions on ^^G and about a factor of 5 
larger for the neutral-current reactions [323] (although they can be induced by all neutrino flavors). 

6.2 Oxygen 

The double-magic nucleus is a potential neutrino target in water Gerenkov detectors like Su- 
perkamiokande. Like ^^G it has very large Q values for charged-current reactions. Furthermore, the 
closed-shell structure of the ground state strongly suppresses GT transitions (in the IPM limit the 
GT strength is zero). These two effects strongly reduce the neutrino cross sections on for supernova 
neutrinos. In fact, charged-current reactions have been estimated to give a negligible event rate in water 
Gerenkov detectors for galactical supernovae with a distance of 10 kpc or 30000 light years (the distance 
to the galactical center) [324, 325]. Due to their expected higher average energies neutrinos might 
excite the spin dipole resonances in which will then decay by emission of protons or neutrons. If 
this decay leaves the daughter nucleus, or ^^ 0 , in particle-bound excited states (i.e. at excitation 
energies E^ < 10.7 MeV in and E^ < 7.3 MeV in ^^O) the subsequent 7 decay cascade of these 
excited states might be observable in Superkamiokande, in particular as both daughter nuclei have no 
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Figure 29: (left) Comparison of cross sections on calculated in shell model and RPA 

approaches (from [321]); (right) Contributions to the (z/e,e“) cross sections on ^°Ar from the IAS and 
from allowed (calculated in shell model) and forbidden (calculated by RPA) transitions (from [328]). 


excited state below 5 MeV. In Ref. [326] this decay scheme has been estimated to yield about 700 events 
in Superkamiokande for a supernova in the galactical center. However, the estimate has been made for 
a Vx Fermi-Dirac neutrino spectrum with temperature T = 8 MeV (and zero chemical potential). More 
modern supernova simulations indicate smaller average energies (temperatures) for z/j, neutrinos, which 
would noticeably reduce the expected event rate. The main observational channel of supernova 
neutrinos in water Cerenkov detectors appears to be coherent elastic scattering on protons [318]. 

As the Gamow-Teller contribution is very small, the charged- and neutral-current cross sections 
on are dominated by spin-dipole transitions. These have recently been modelled in a shell model 
calculation [321] considering the p and sd shells and modifying the residual interaction to account for 
effects of the tensor force [327]. This calculation hnds the major dipole strength for the 2“ channel with 
excitations around = 15 and 22 MeV. As is demonstrated in the left panel of Fig. 29, the calculated 
shell model ^®0(z/e, e“)^®F cross sections agrees very well with results obtained within the Continuum 
RPA [190] for neutrino energies up to 100 MeV. 

6.3 Argon 

Liquid argon detectors are large time-projection chambers with a very good energy resolution and full 
particle track reconstruction [191]. For supernova neutrinos the main goal is to observe Ve neutrinos 
from ^°Ar(z/e, e“)‘^°K* by tagging on the gamma’s from the deexcitation of the hnal states in [191]. 
A similar detection scheme is possible for the observation of Ve neutrinos. However, as we will mention 
below, the respective cross section is reduced due to the suppression of Gamow-Teller distributions. 
The ICARUS detector is a currently operational liquid argon supernova neutrino detector [329, 330]. 

The ^°Ar(z/e, e“)^°K cross section for supernova neutrinos is dominated by the Fermi contribution 
to the IAS in at = 4.38 MeV and the Gamow-Teller transitions. The hrst can readily be 
calculated applying the {N — Z) sum rule for Fermi transitions, while the second is constrained by 
the recent measurement of the B(GT_) strength of ^“^Ar by a (p, u) charge-exchange experiment [331]. 
However, it has been found that the GT strengths assigned to transitions to two low-lying states at 
excitation energies E^ = 2.289 MeV and 2.730 MeV is relatively large, hence signihcantly contributing 
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to the charged-current cross section, but it is noticably different than the B(GT) value derived from 
the decay of the isospin-analog ^°Ti [332], This difference has been recently investigated in a (^He,t) 
experiment, which has determined the relative strength of the two transitions and found a value close 
to the (p, n) experiment and in clear contrast to the /d-decay data (quoted in [333]). A shell model 
calculation performed in the pf-sd model space has nicely reproduced the B(GT) strength determined 
in the {p,n) experiment, including the relative transition strengths to the two low-lying states [328]. 
Recently a shell model calculation, including a residual interaction with consideration of Goulomb 
contributions, hence breaking isospin symmetry, has explored possible nuclear structure reasons for the 
differences between the B(GT) values derived from the {p,n) measurement and the '^^Ti fd decay [333]. 

The authors of Ref. [333] recommend to use the Gamow-Teller strength determined from the (p, n) 
experiment to derive the '^°Ar(i/e, e“)^°K cross section for supernova neutrinos. This has indirectly been 
done by Suzuki and Honma, as their shell model Gamow-Teller strength distribution gives a very good 
account of the one derived from the (p, n) data [328]. The calculation shows that there is additional 
GT strength outside of the experimental energy window explored by the (p, n) data. While theory and 
experiment give a B(GT) strength of about 5 units up to excitation energies of 8 MeV in the 
calculated total strength is 7.21 units. The recent shell model calculation also gives a better account 
of the data than obtained in a model space restricted to 2p-excitations from the sd to the pf shell 
[334]. Based on a hybrid model approach (shell model GT strength distribution, the transition to the 
IAS from the Fermi sum rule and forbidden transition calculated within the RPA for multipolarities 
A = 1-5 Suzuki and Honma have calculated the ‘^°Ar(z/e, cross section for neutrino energies up 

to El, = 100 MeV (see right panel of Fig. 29). As expected, the cross section is dominated by allowed 
contributions up to Ei^ = 50 MeV, which covers the range of supernova Ue energies. It is interesting to 
note that the cross sections obtained within the hybrid model agrees well with calculations where all 
multipole contributions have been derived within the RPA [18] except at low neutrino energies where 
the hybrid model cross sections are larger by about 20-40 %. We note again that at low neutrino 
energies the cross section is sensitive to the detailed GT strength distribution which, in agreement with 
recent data, is well described by the shell model. Hence the results presented in Ref. [128] should be 
more reliable than the RPA results of ref. [18]. The neutrino-induced charged-current reactions have 
also been studied within the Local Density Approximation taking account of Pauli blocking and Fermi 
motion [335]. The results agree reasonably well with those from the hybrid model and the RPA approach 
for neutrino energies Ei, > 20 MeV. Gross sections presented on the basis of a QRPA calculation [336] 
deviate from the values obtained in the hybrid model and the RPA study by about a factor of 3. 

In the IPM, the GT+ strength for ^‘^Ar (with 18 protons and 22 neutrons) vanishes as all GT 
transitions, in which a proton is changed into a neutron, are Pauli blocked. Gross-shell correlations might 
introduce a non-vanishing GT+ strength, but it is expected to be small. Based on this assumption, the 
^°Ar(z/e, e+)^°Gl cross section has been calculated solely within an RPA approach [18]. Due to the missing 
or highly suppressed allowed contributions, the (Pe, e’*') cross section on “^^Ar is much smaller than the 
cross section for supernova neutrino energies [18], making the observation of Pg neutrinos by 
liquid argon detectors a challenge. It would be still desirable to put the low energy (Pg, e^) cross section 
on a more reliable basis by calculating the allowed and forbidden transitions to low-lying states within 
shell model calculations which take cross-shell correlations properly into account. Such studies have 
been proven to be up to the task in studying the Ml strength in ^®Ar [171] or the isotope shift in calcium 
[170]. 

6.4 Lead 

The Fermi and Ikeda sum rules both scale with neutron excess (N — Z). As the charged-current response 
induced by supernova z/g neutrinos (with average energies around 12 MeV) are dominated by Fermi and 
GT transitions, the charged-current cross sections on lead is expected to be significantly larger than 


on other materials like iron or carbon. This makes attractive target for snpernova nentrino 

detectors like the HALO detector in Canada [337]. The main detection schemes are the (z/e,e“) and 
{ux, h'x) reactions where in both cases the final nnclei, ^°®Bi and ^°®Pb, respectively are left in excited 
states which mainly decay by particle emission. Dne to the relatively low excitation energies involved 
and the large Coulomb barriers associated in the proton channel, the particle decay occurs mainly by 
emission of (one or two) neutrons. In HALO these neutrons can be detected from capture on ^He, after 
they have been moderated in polypropylene [337, 191]. Hence the calculation of the partial neutrino- 
induced neutron spallation cross sections on ^°^Pb are of relevance to interpret supernova neutrinos 
detected by HALO. 

There have been several evaluations of the (z/e,e“) cross sections on ^‘^®Pb, which is the dominating 
isotope of natural lead. Unfortunately the calculation of the allowed contributions to the cross sections 
within the shell model is yet computationally not feasible. The Fermi contribution to the cross section 
can be derived from the transition to the known IAS in ^°®Bi and the Fermi sum rule B{F) = {N — Z). 
All studies have calculated the contributions from the other multipoles on the basis of RPA variants 
which reproduce the strength and the position (around = 15 MeV in ^°®Bi) of the Gamow-Teller 
distribution as well as of the spin-dipole resonances [165, 338, 339] The B(GT_) strength in ^ospp jg 
hxed by the Ikeda sum rule, however, it has been renormalized in Refs. [165, 339] (by (0.7)^ and (0.8)^, 
respectively) to account for the quenching of the GT strength at low excitation energies, while Ref. 
[338] used the unquenched qa value. Despite these differences the calculations give quite similar cross 
sections (within 20 — 30%) as function of neutrino energy. The results are, however, noticeable smaller 
(by a factor 3) than those obtained in the first study of the reaction on ^ospp [ 340 ] which is 

due to the differences in the nuclear model, the treatment of the multipole operator expansion and 
the Goulomb correction for the electron in the outgoing channel [165, 338]. A QRPA calculation, on 
top of a Hartree-Fock study with BGS pairing, has been presented in Ref. [160]. The obtained results 
are about a factor of 2 larger than in Refs. [165, 339] (also for ^®Fe), likely caused by the neglect 
of renormalization of the GT contribution to the cross section. Suzuki and Sagawa presented (z/g, e~) 
cross sections which have been obtained using GT data from a (p, n) experiment performed at Osaka 
[341] and adjusting their Hartree-Fock -|- Tamm-Dancoff approach for the hrst-forbidden response to 
the peaks of the spin-dipole resonances [342]. Furthermore, the spreading and quenching of the GT 
response has been considered by coupling to 2p-2h conhgurations. Although the calculations have been 
performed assuming a muon-decay-at-rest rather than a supernova neutrino spectrum, and are thus 
not directly relevant to the present discussion, the result, however, is satisfyingly: Suzuki and Sagawa 
obtain (3.2 x 10“^® cm^) for the cross sections [342], in close agreement to the RPA result of 

Ref. [165] (3.62 x IQ-^^ cm^). 

Jachowicz and co-workers [343, 344] have investigated the influence of different neutrino spectra and 
of neutrino oscillations on the neutrino signal in lead-based supernova detectors. Total charged-current 
cross sections for the chain of lead isotopes for supernova neutrinos, assuming a Fermi-Dirac 
distribution with temperature T = 4 MeV and zero chemical potential, are reported in Refs. [104, 345] 
based on RPA calculations yielding quite similar cross sections. RPA-based studies have been performed 
for the neutral-current {u, E) cross sections on [165, 339]. Again the results agree quite well (within 
10-15 %). 

It might not be surprising that the RPA sudies of Refs. [165, 338, 339] predict similar cross sections, 
as they are constrained by similar indirect data like the positions and strengths of the leading responses. 
However, the calculations differ noticeably in their treatment of the quenching of the GT strength. 
Here a recent measurement of the Gamow-Teller and spin-dipole responses by a (p, n) charge-exchange 
reaction with polarized protons and neutrons has determined the GT_ strength up to excitation energies 
of nearly 50 MeV. The data yield a GT_ strength of about 85 units up to an excitation energy of 25 
MeV, which includes the GT resonance at about 15 MeV. This value is significantly smaller than the 
Ikeda sum rule, 3{N — Z) = 132, corresponding to a quenching of the low-lying GT_ strength by a factor 
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of (0.8)^ as assumed in Ref. [339]. Ref. [104] used too strong a quenching, hence the cross sections might 
be slightly too low, while Ref. [338] did not consider any quenching and might have overestimated the 
GT contribution to the cross section. The GT and spin-dipole data show another known shortcoming 
of RPA calculations as they often do not describe the observed fragmentation and spreading of the 
resonant strength well due to missing correlations. As we have discussed above this is not too relevant 
for the total cross sections if the neutrino energies are sufficiently large. However, these shortcomings 
are quite relevant if one wants to describe the neutron spallation cross sections as the GT resonance 
resides quite closely to the threshold of the 2n decay channel which opens at = 14.98 MeV. Such 
studies have been performed for example in Refs. [165, 339] where the latter simply assumes that all 
states with excitation energies between 100 keV above the In threshold at = 6.9 MeV to 2.2 MeV 
above the 2 n threshold decay by emission of one neutron, while states at higher excitation energies 
decay by multiple neutron emission. Ref. [165] treats the various decay modes of excited states in ^°®Bi 
as a cascade of decays followed on the basis of the statistical model. Ref. [165] uses the information from 
the statistical decays to calculate neutron spectra for various supernova neutrino spectra. These results, 
however, should be effected by the underestimation of the spreading and fragmentation of the resonant 
responses (see the discussions in Refs. [165, 18]). These references discuss whether the observation of the 
neutron count rate in a lead supernova detector might be a suitable tool to study neutrino oscillations 
as proposed in [340]. 

6.5 Molybdenium and cadmium 

There have been recently several studies of neutral- and charged-current reactions on molybdenium and 
cadmium isotopes induced by supernova neutrinos. These investigations have been motivated to explore 
the sensitivity of the Gadmium Zinc Telluride 0-Neutrino Double-Beta Research Apparatus (GOBRA) 
experiment [346] and the Molybdenium Observatory of Neutrinos (MOON) [347] to the observation 
of supernova neutrinos. The GOBRA and MOON detectors are both designed and planned for the 
detection of neutrino-less double beta decay. If this mode is successfully observed it will establish the 
Majorana nature of the neutrino and, even more importantly, will establish violation of lepton number 
conservation (e.g. [348]. The GOBRA detector is made out of GdZnTe crystals which serve as both 
source and detection material. Its potential to detect neutrino-less double beta decay rests mainly on 
the two nuclides ^^®Gd and ^^‘^Te. The MOON detector is a tracking-calorimeter detector which is highly 
enriched in ^‘^^Mo [349]. Besides searching for neutrino-less double beta decay the MOON detector is 
also anticipated as solar neutrino detector. 

Motivated by the relevance for double-beta decay, GT strength distributions for ^^‘^Mo and ^^®Gd 
have been determined by charge-exchange reactions. Using a (^He,t) experiment at RGNP Thies et al 
showed that the GT+ strength on ^‘^^Mo is basically concentrated in a single transition to the ^°°Nb 
ground state [350], likely corresponding to the transformation of a gg/ 2 - proton into a 5 ^ 7/2 neutron. The 
GT strength between these two ground states has independently been deduced from the ^‘^‘^Tc electron 
capture branching ratios, obtained at the IGISOL facility [351]. We note that a similar domination of 
the GT strength by a single transition (to the state at E^ = 0.69 MeV) has been observed in the 
GT+ strength of ®®Mo, deduced by (d,^ He) and (^He,f) charge-exchange experiments. For ^^®Gd the 
GT_ strength has been determined recently using the {p,n) reaction [352, 353]. The centroid of the 
strength associated with the GT resonance is located at an excitation energy around E^ = 14 MeV, 
with some additional strength forming a shoulder in the GT distribution around E^ = 8 MeV. The 
GT_ strength on ^^®Gd had previously been determined by a (^He,t) experiment [354] which turned 
out, however, to be erraneous due to the use of a natural cadmium target [355]. The GT strength 
between the ^^®Gd and ^^®In ground states has also been determined from the electron capture decay 
branch of ^^®In [356] yielding a transition strength in agreement with the value deduced from the {p, n) 
charge-exchange experiment [353]. 
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The experimental GT strength distributions are valuable constraints for calculations discussing the 
MOON and COBRA detectors as devices to observe supernova neutrino. A detailed study of neutrino- 
and antineutrino scattering off ^^®Cd has recently been reported in Ref. [357] performed within the 
charge-changing mode of the QRPA. The authors performed detailed studies of the GT_ strength dis¬ 
tribution for several Skyrme interactions supplemented by hnite-range (isovector and isoscalar) pairing 
interactions and compare it to the results obtained with the Bonn one-boson exchange potential. (Sim¬ 
ilar studies based on the Bonn potential have been previously reported in [358]. These authors adjusted 
the parameters of the force for each multipole individually.) Only the Bonn potential was found to 
reproduce the giant GT resonance at 14 MeV and its satellite around 8 MeV, while the calculations 
performed with the Skyrme forces missed these energies typically by 1-2 MeV [357]. These differences in 
the position of the strong GT transitions translate into different predictions for the {Uf,, e~) cross section 
on ^^®Cd for supernova neutrinos, which, as calculated in Ref. [357], is dominated by GT transitions. 
Using the cross section obtained with the Bonn potential as benchmark, the Skyrme interactions yield 
results which can deviate upto a factor of 2; however, most Skyrme interactions predict cross sections 
which agree within 50%. For the (z/g, e"*") cross section on ^^®Gd for supernova antineutrinos the cal¬ 
culations predict results which agree within a factor of 2 among each other, with the Skyrme forces 
typically yielding a somewhat larger value than the Bonn interaction. The largest contribution to the 
cross sections arizes from hrst-forbidden transitions to the J = 1~ states in ^^®Ag. The partial cross 
sections reported in Ref. [358] for the (differently parametrized) Bonn potential are very similar, except 
for the contribution of the J = O’*" multipole to the cross section which is noticeably larger 

in the study of [358] than found in Ref. [357], in particular compared to the results for the Skyrme 
interactions. 

The formalism of charge-changing QRPA has also been adopted to calculate the cross sections for 
neutral-current reactions on stable cadmium and molybdenium isotopes induced by supernova neutrinos 
of all flavor. For the cadmium isotopes the inelastic cross sections for supernova neutrinos are 
dominated by the contributions from the J = multipole (GTq strength) [359]. For the odd-A 
isotopes Almosly et al. hnd a surprisingly large vector contribution to excited J = O’*" states 

[359]. 

Also the calculation of neutral-current reactions on the stable molybdenium isotopes show a some¬ 
what larger cross sections for the odd-A nuclei ®®’®^Mo compared to their even-even neighbor iso¬ 
topes [360]. The difference amounts to about 40% for neutrinos with energies typical for supernova 
z/g’s and decreases with neutrino energy. The neutral-current response of supernova neutrinos on odd-A 
and even-even molybdenium isotopes has also been studied in Refs. [361, 362]. Interestingly the au¬ 
thors hnd a relatively large contribution of vector transitions to excited O’*" states for the odd isotopes 
95,97 mo [361], which is less pronounced for the even-even isotopes [362]. The charged-current reactions 
on the odd isotopes ®®’®^Mo for supernova neutrinos has been explored in Ref. [363]. The authors 
hnd the (z/g,e“) cross section dominated by allowed Fermi (to the IAS) and GT transitions, while the 
(z/g, e+) cross section has some forbidden contributions arizing from 1“ and 2“ multipoles, besides the 
dominating GT contribution. 

Engel et al. [364] have studied the charged-current cross sections on ^°°Mo, replacing in the domi¬ 
nating 1+ channel the QRPA results with the GT strength measured by the (^He, t) charge-exchange 
reaction [354]. These authors discuss how the MOON detector might distinguish between charged- and 
neutral-current events, assuming that neutral-current reactions will excite the molybdenium isotopes to 
states which mainly decay by emitting neutrons and successive y-rays [364]. Hence it appears useful if 
the studies reported in Refs. [361, 362] would be extended to calculations of the partial cross sections 
with particle emission in the hnal channel. 
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7 Summary 


Neutrinos have been identified as one of the essential players in core-collapse supernovae. This insight 
- predicted by models - has been unambiguously verified by the observations from supernova SN1987A. 
The undisputably most important role of neutrinos for the supernova dynamics concerns the energy 
balance: due to their tiny cross sections with the surrounding matter neutrino emission efficiently cools 
the inner core of the collapsing core and, after bounce, carries away an overwhelmingly large portion 
of the gravitational binding energy released in the explosion. Even despite the tiny cross sections, the 
tremendously large number of neutrinos involved makes them an important means of energy transport 
during the explosion, reviving the stalled shock front. 

Neutrinos also play a crucial role for the explosive supernova nucleosynthesis as their interaction 
with the hot matter surrounding the protoneutron star (mainly free nucleons, but also light nuclear 
clusters) determines the ratio of protons and neutrons available for synthesis of heavier nuclei in the up 
process or potentially the r-process. While (anti)neutrino absorption on protons is at the very heart 
of the up process, the role of neutrinos during the mass flow to heavy nuclides in the astrophysical r- 
process (if it occurs in the core-collapse supernova environment) is likely to be small. However, particle 
spallation induced by neutrino reactions on abundant nuclei in the outer layers of the star has been 
identified as the synthesis mechanism of selected isotopes like ^^®La and ^®°Ta in a process called 
neutrino nucleosynthesis. 

In several of these aspects of supernova dynamics and nucleosynthesis, neutrinos interact with free 
nucleons. The respective cross sections are well defined based on the framework of the electroweak 
theory. However, in the dense and hot supernova environment these interactions are modified by 
correlations which is a topic of ongoing research interest. Neutrino interactions with nuclei can occur 
during collapse as well as after bounce during the various nucleosynthesis processes. During collapse, 
an important neutrino process is coherent scattering on nuclei which leads to neutrino trapping in the 
final phase of the collapse. Except for small corrections induced by isospin breaking and additional 
Gamow-Teller contributions (at finite temperature and for odd-H and odd-odd nuclei), the coherent 
cross sections are determined by the Fermi transition to the Isobaric Analog State and hence well under 
control. Other neutrino processes during collapse include absorption on nuclei and inelastic neutrino 
scattering with nuclei. As electron capture on nuclei is the other very relevant weak-interaction process 
during collapse, mainly u^ neutrinos are present in the core. They can be absorbed on neutrons in nuclei. 
As this is the inverse reaction of electron capture, the respective cross section is derived applying the 
principle of detailed balance. Hence, its accuracy has benefitted from the improved description of 
electron capture on nuclei made possible in recent years by more precise data from advanced charge- 
exchange experiments and from progress in nuclear modelling, mainly due to the application of shell 
model techniques. 

The shell model has also been the many-body technique which has been applied in the first derivation 
of inelastic neutrino-nucleus scattering cross sections at finite supernova temperatures. The calcula¬ 
tions have been validated against precision data derived from inelastic electron scattering. As these are 
taken for nuclear ground states only, the validation is restricted to cross sections at temperature T = 0. 
However, neutrino-induced deexcitation of nuclear states, thermally populated at finite temperatures, 
are found to strongly increase the inelastic cross sections at low neutrino energies. Although the shell 
model is expected to describe also the (Gamow-Teller) transitions on excited states adequately well, a 
full state-by-state summation to obtain the inelastic cross section at supernova temperatures is compu¬ 
tationally unfeasible and some approximation has to be applied. However, it has been quite gratifying 
that calculations performed within the Thermal QRPA approach confirmed the strong increase of the 
cross sections at low neutrino energies at finite temperature. This method derives the cross sections 
consistently from a thermal ensemble and accounts for effects induced by the melting of nucleon pairs 
at finite temperature. It is important that both approaches - the TQRPA and the shell model - predict 
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very similar cross sections at higher neutrino energies as supernova simulations have indicated that in¬ 
elastic scattering of such high-energy neutrinos off nuclei leads to down-scattering of the neutrinos and 
hence strongly reduces the tail of the neutrino spectra produced in the neutrino-burst just after bounce. 
As a consequence, this change in spectra reduces the expected event rates for the future observation of 
neutrinos from the burst phase by supernova neutrino detectors. Inelastic neutrino-nucleus scattering 
is also a mean of energy exchange between neutrinos and matter during collapse. But simulations have 
indicated that this process only adds rather mildly to neutrino thermalization, which is dominated by 
inelastic neutrino scattering on electrons. 

Neutrino types other than electron neutrinos are produced during infall by deexcitation of thermally 
populated nuclear states by neutrino-pair emission. However, for the supernova dynamics this process 
has been recently shown as unimportant as the cross section for z/g production by electron capture 
is several orders of magnitude larger and dominates the energy balance of the core during collapse. 
This finding is based on a reasonable, but quite schematic description of the neutrino-pair deexcitation 
process, but is unlikely to change if more accurate cross sections, for example derived on the basis of 
the Thermal QRPA model, become available. 

Neutrino-induced reactions have been proposed to play roles in various supernova nucleosynthesis 
processes; e.g. in the up and r-process nucleosynthesis and obviously in the process called neutrino 
nucleosynthesis. To allow for systematic and detailed investigations into the importance of neutrino- 
nucleus reactions for these various processes, a global set of partial differential cross sections for neutrino- 
induced charged- and neutral-current reactions on nuclei (up to charge number Z = 78) has been 
derived. Based on the observation that for the moderate supernova neutrino energies the initial neutrino- 
nucleus excitation function is dominated by the contributions of giant resonances, the calculation of 
the respective cross sections has been performed within the Random Phase Approximation, which is 
constructed to give a fair account of the energy and the total transition strength of these collective 
resonances. In these initial processes the neutrino excites the nucleus to levels which subsequently 
decay by particle or 7 emission. The decay of the excited nuclear levels is followed within the Statistical 
Model considering a cascade of multi-particle decays. 

Inserting these neutrino-nucleus cross sections into supernova nucleosynthesis simulations indicates 
that neutrino reactions on nuclei have only limited impact on the calculated abundances. For the up 
process Ue absorption on protons bound in nuclei can in principle compete with the absorption on free 
protons, but is kinematically suppressed due to the relatively large Q values of the neutron-deficient 
nuclei encountered in the up process mass flow. Current supernova simulations indicate only a ‘weak 
r-process’ which produces nuclides up to the second r-process peak (mass numbers A < 130). In 
this scenario neutrino-nucleus reactions are found unimportant. Uncertainties related to supernovae 
as r-process sites are introduced due to the yet insufficiently understood roles of collective neutrino 
oscillations and of sterile neutrinos as well as to the undetermined mass hierarchy of neutrinos. 

More than two decades ago, neutrino-induced spallation of abundant nuclei in the outer layers of 
a supernova has been identified as a major nucleosynthesis source of selected nuclides. This finding 
has been confirmed over the years using improved stellar models, more accurate neutrino-nucleus cross 
sections and, last, but not least, refined spectra for the various types of supernova neutrinos. Although 
the predictions for these spectra shifted systematically to lower average energies (or temperatures if 
one describes the spectra as usual by a Fermi-Dirac distribution with zero chemical potential) stellar 
models still imply that the nuclides ^Li, ^^B, ^^®La and ^®°Ta are produced by neutrino nucleosynthesis 
nearly in solar abundances. As ^^B is the product of neutral-current reactions on induced by 
and u-r neutrinos and their anti-particles and ^^®La and ^®Ta are made by (z/e,e“) reactions on ^^®Ba 
and respectively, neutrino nucleosynthesis can serve as a thermometer for the spectra of those 

neutrino types, which have not been observed from supernova SN1987A. 

The observation of neutrinos from a future close-by supernova, preferably from our galaxy, is the 
ultimate goal of the various neutrino detectors either operational or under construction. If the supernova 
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is indeed close enough, e.g. from the galactic center, these detectors have the ability to test the hierarchy 
of neutrino energies predicted by supernova models. Furthermore they will have the time resolution 
to probe the evolution of neutrino luminosities during collapse and explosion, including the prominent 
burst of electron neutrinos expected just after bounce when heavy nuclei are dissociated into free 
nucleons and the electrons, which have survived the collapse, are captured on free protons. To translate 
the observed neutrino event rates from a supernova requires the knowledge of the respective neutrino- 
induced reactions on the detector nuclei. These cross sections have been derived on the basis of nuclear 
models currently available to treat the related many-body problem. Hence for lead or cadmium targets 
only cross sections calculated within the frame work of the RPA are available, while for lighter target 
nuclei more sophisticated nuclear models have been applied. As a consequence for the lighter nuclei 
a more accurate description of the detailed nuclear response has been considered than for the heavier 
nuclei, which is of relevance if particle decays of levels in the daughter nucleus might serve as detection 
scheme. Furthermore some more fundamental questions like the quenching of higher multipoles, in 
relation to the one established for shell model calculations of the Gamow-Teller strength when compared 
to data, or the role of two-body currents has to be resolved. 

Even if a few questions concerning the role of neutrino-induced reactions on nuclei for the supernova 
dynamics and nucleosynthesis remain insufficiently answered, the held has witnessed signihcant progress 
in recent years. The decisive push in this held in specihc, and, of course, to our general understanding 
of supernovae globally would come from the observation of a nearby supernova by 7 ’s, neutrinos and 
the composition of the ejecta. Perhaps Nature is so kind and we do not have to wait for too long. The 
observers are ready. 
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